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ABSTRACT 

We consider the interaction of a relativistically-moving shell, composed of thermal photons, a reversing mag- 
netic field and a small admixture of charged particles, with a dense Wolf-Rayet wind. A thin outer layer of 
Wolf-Rayet material is entrained by the jet head; it cools and becomes Rayleigh-Taylor unstable, thereby provid- 
ing an additional source of inertia and variability. The gamma-rays emitted by the shell load the ambient medium 
with electron-positron pairs and, close to the engine, force this medium to move at nearly the same speed as the 
shell. We show that this pre-acceleration of the wind material defines a characteristic radiative compactness at the 
point where the reverse shock has completed its passage back through the shell. We argue that the prompt gamma- 
ray emission is triggered by this external braking, at an optical depth ~ 1 to electron scattering. Torsional waves, 
excited by the forced reconnection of the reversing magnetic field, carry a fluctuating current, and are damped 
at high frequencies by the electrostatic acceleration of electrons and positrons. We show that inverse Compton 
radiation by the accelerated charges is stronger than their synchrotron emission, and is beamed along the magnetic 
field. Thermal radiation that is advected out from the base of the jet cools the particles. The observed relation 
between peak energy and isotropic luminosity - both its amplitude and scaling - is reproduced if the blackbody 
seeds are generated in a relativistic jet core that is subject to Kelvin-Helmholtz instabilities with the Wolf-Rayet 
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envelope. This relation is predicted to soften to £p e ak ~ L^ below an isotropic luminosity L[ so 
Spectrally harder bursts will arise in outflows which encounter no dense stellar envelope. The duration of spikes 
in the inverse-Compton emission is narrower at higher frequencies, in agreement with the observed relation. The 
transition from prompt gamma-ray emission to afterglow can be explained by the termination of the thermal X-ray 
seed and the onset of synchrotron-self-Compton emission. GLAST will probe the mechanism of particle heating 
at the reverse shock by measuring the inverse-Compton scattering of the seed photons. 

Subject headings: gamma rays: bursts; supernovae: general; turbulence; radiation mechanisms: non-thermal 
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1. INTRODUCTION 

Models of gamma-ray burst (GRB) outflows have long been 
plagued by uncertainties in the mechanism by which gamma- 
rays are emitted, including the underlying instability that trig- 
gers the rapid variability, the mechanism by which particles are 
energized, and the location of the gamma-ray emitting region. 
These uncertainties are closely interlocked. Because there is 
no shortage of plausible mechanisms which could, in principle, 
generate gamma-rays in a relativistic outflow, the lack of con- 
vergence on one issue has tended to impede progress on all the 
others. 

The data offer some important clues which have not yet been 
fully exploited by theoretical models. For example, individ- 
ual GRB light curves show only weak changes in the pattern 
of variability from the beginning to the end of the burst (Be- 
loborodov, Stern, & Svensson 1998; Ramirez-Ruiz & Fenimore 
2000). This, in spite of the basic requirement - based on simple 
considerations of gamma-ray opacity (Cavallo & Rees 1978) - 
that the non-thermal gamma-rays must be emitted in a region 
that is vastly larger than the central engine - by some 6-8 or- 
ders of magnitude. One deduces that this dissipative zone must 
be concentrated within a fairly narrow range of radius, in spite 
of the large distance that the outflow must travel. In addition, a 
well-defined correlation has emerged between the peak energy 
of gamma-ray bursts and their isotropic luminosity, in a sam- 
ple of BeppoS AX and HETE-II bursts with measured redshifts 
(Amati et al. 2002; Ghirlanda et al. 2004; Lamb, Donaghy, & 
Graziani 2005). A third clue is that the the short gamma-ray 
bursts have systematically harder spectra than the long bursts. 



The recent burst GRB 050509b had and isotropic luminosity of 
~ 3 x 10 49 ergs s" 1 at a redshift of z = 0.2, but a much harder 
spectrum than the Amati relation would imply (Bloom et al. 
2005). This observation has interesting implications for the lo- 
cation of the dissipative zones in a GRB outflow. 

Strong circumstantial evidence has emerged that long GRBs 
involve the ejection of relativistically moving material from ~ 
2-3M Q black holes, which are formed in the collapse of mas- 
sive stars (Woosley 1993; Paczyhski 1998). The gamma-ray 
burst 030329 coincided with the energetic Type Ic supernova 
2003dh (Stanek et al. 2003). The optical spectrum of this su- 
pernova was similar to that of the energetic Type Ic SN 1998bw 
(Galama et al. 1998). This previous supernova coincided with a 
gamma-ray burst (GRB 980425) of a much smaller isotropic lu- 
minosity, which may therefore have involved a buried jet. The 
progenitors of the long GRBs are apparently Wolf-Rayet stars 
surrounded by fast and dense winds. The origin of the short 
gamma-ray bursts is more uncertain, and may involve a mix 
of sources including binary neutron star systems (Eichler et al. 
1989) and extragalactic Soft Gamma Repeaters (Duncan 2001). 

The Wolf-Rayet/jet model offers two specific locations 
where the physics of gamma-ray emission can be well con- 
strained: first, where the collimated relativistic outflow breaks 
out of the Wolf-Rayet star; and, second, where it begins to 
be decelerated substantially by its interaction with the external 
medium (including the outermost thin layer of stellar material 
that is entrained during the breakout of the relativistic jet). Our 
first goal in this paper is to define how the deceleration phase is 
modified by pair creation in the Wolf-Rayet wind, and the as- 
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sociated radiative acceleration of the pair-loaded wind material. 
Previous treatments of this process (Thompson & Madau 2000; 
Meszaros, Ramirez-Ruiz, & Rees 2001; Beloborodov 2002) as- 
sume a given pattern of gamma-ray emission, and do not exam- 
ine its feedback on the forward and reverse shocks at the very 
earliest stage where the prompt burst is generated. One finds 
that the relativistic ejecta have a characteristic compactness at 
the point where the reverse shock wave completes its passage 
through the ejecta shell, because the process of pre-acceleration 
shuts off below a critical value of the fiuence of the gamma-rays 
flowing across the forward shock. Rapid deceleration of the rel- 
ativistic ejecta follows once this threshold is reached. The com- 
pactness is, in fact, regulated to a value somewhat larger than 
unity, and the optical depth through the swept-up e + -e~ pairs 
has a similar value. One therefore obtains a direct feedback be- 
tween the physical conditions in the dissipative zone, and the 
radiative mechanism. 

The gamma-ray emission depends on the presence of some 
irregularities in the outflow. These irregularities could be 
caused by stochastic flips in the sign of a dynamo-generated 
magnetic field that is swept into the outflow (Thompson 1994); 
or, alternatively, they could involve fluctuations in momen- 
tum and energy flux that manifest themselves as internal shock 
waves in a relativistic particle outflow (Rees & Meszaros 1994). 

The seed irregularities are, in both of these models, im- 
printed into the outflow near its base. Nonetheless, the mod- 
els are distinguished by the isotropy of the emission in the 
bulk frame. The synchrotron or inverse-Compton radiation of 
shock-accelerated particles is nearly isotropic in the rest frame 
of the shocked fluid. The observation of fast variability in many 
gamma-ray bursts then requires that the gamma-rays be emit- 
ted far inside the radius at which the bulk motion begins to 
be reduced significantly (Sari & Piran 1997). The basic dis- 
advantage of this approach is that the location of the dissipative 
zone depends sensitively on the frequency and amplitude of the 
seed irregularities. Calculating the pattern of these irregulari- 
ties from first principals requires a detailed theoretical under- 
standing of almost every aspect of the central engine and its 
outflow. In spite of recent herculean efforts to simulate jet out- 
flows (McKinney & Gammie 2004; de Villiers, Staff, & Ouyed 
2005; McKinney 2005a,b), the model retains considerable free- 
dom in practice. 

A hybrid model of jet dissipation is suggested by consider- 
ing the influence of jet breakout on the subsequent motion of the 
reverse shock wave. A thin outer layer of the Wolf-Rayet star 
accumulates energy from the relativistic beam emerging from 
below, and is pushed outward to higher speeds. Existing calcu- 
lations of this 'breakout shell' are limited to a one-dimensional 
approximation (Waxman & Meszaros 2003). We calculate the 
column density of this shell, and the growth of its Lorentz fac- 
tor, by taking into account the slippage of the swept-up matter 
to the side of the beam. Rayleigh-Taylor instabilities allow the 
relativistic fluid to emerge through the breakout shell in isolated 
spots, which can radiate independently. 

Our second goal in this paper is to establish a realistic mech- 
anism by which the gamma-ray emission may be triggered 
directly by the deceleration off the 'breakout' shell and the 
external wind medium. Calculations of prompt high-energy 
emission from the forward shock (Meszaros & Rees 1993; 
Meszaros, Rees, & Papathanassiou 1994; Chiang & Dermer 
1999; Sari & Esin 2001) have neglected the effects of pre- 
acceleration of the ambient medium, have assumed that the 



radiating particles are shock-accelerated, and therefore cannot 
easily account for fast variability in the emission. This forces 
us to focus on an alternative source of energy in the outflow - 
in particular, a strong non-radial magnetic field, which is inde- 
pendently motivated by considerations of the launching of the 
relativistic outflow. The magnetic field provides a character- 
istic direction for particle acceleration, and provides plausible 
instabilities which cause impulsive and localized injections of 
energy. 

The relative proportions of energy carried by Poynting flux 
and by ions are still very uncertain in observed extragalactic jet 
sources (e.g. Sikora et al. 2005). Basic conservation laws limit 
the transfer of energy from the magnetic field to the particles 
through internal instabilities, if the inner boundary conditions 
of the outflow are time-invariant. Magnetic flux, in particu- 
lar, is not easily eliminated from a fast-moving astrophysical 
fluid with a very high magnetic Reynolds number. The ion lu- 
minosity can increase by a logarithmic factor with respect to 
the Poynting luminosity in a smooth outflow (Begelman & Li 
1994). The magnetic field will, nonetheless, have a strong soft- 
ening effect on internal shock waves if the ordered Poynting 
luminosity Lp comprises 10 percent of the outflow luminosity 
(Kennel & Coroniti 1984; Zhang & Kobayashi 2005); and will 
contribute substantially to the post-shock pressure if the Poynt- 
ing luminosity is somewhat smaller than this. Recent jet sim- 
ulations suggest an approximate equipartition between kinetic 
and magnetic energy at a distance of ~ 10 10 cm from a stellar- 
mass black hole (McKinney 2005a,b). 

The particle luminosity could also be increased by recon- 
nection of a reversing magnetic field. This process is much 
more effective when the field is being pushed into the external 
medium, at a large distance from the engine where the sound- 
crossing time of the shell is smaller than the flow time (Thomp- 
son 1994; see Lyubarsky & Kirk 2001 for a more detailed dis- 
cussion in the context of pulsar wind termination shocks). Tan- 
gling of the magnetic field allows a reverse shock wave to form 
even when the fast mode speed is higher than the speed of the 
outflow relative to the contact discontinuity. 

Radiation that is advected out from the base of the outflow 
has a characteristic energy of <~ 1 MeV (Goodman 1986). Seed 
black body photons will be present in combination with a mag- 
netic field in many models of the central engine, and therefore 
can play an important role as a source of inverse-Compton seeds 
at a large radius, and help to define the spectral peak of the 
gamma-ray emission (Thompson 1994). Our third goal in this 
paper is to investigate this effect in more detail, in light of more 
recent ideas that the propagation of a relativistic jet through the 
envelope of a Wolf-Rayet star (Zhang, MacFadyen, & Woosley 
2003; Matzner 2003) can create a lower-energy thermal com- 
ponent in the jet cocoon (Ramirez-Ruiz, MacFadyen, & Laz- 
zati 2002). We find that the observed relation between peak 
energy and isotropic luminosity in BeppoSAX and HETE-II 
bursts (Amati et al. 2002; Lamb et al. 2005) arises directly - 
without free parameters - from dissipation in a relativistically- 
moving jet core that is heated but not fully decelerated by 
Kelvin-Helmholtz instabilities with the jet cocoon. (The co- 
coon acquires a similar temperature, but its specific entropy is 
not realistically high enough to seed the gamma-ray emission 
proper.) Instead the portion of the jet material that acquires a 
high entropy within the Wolf-Rayet star is preferentially accel- 
erated to high Lorentz factors after breaking out of the stellar 
photosphere. One can therefore establish a direct connection 
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between the presence of a fireball component at the base of the 
jet, and the observability of that portion of the jet as a gamma- 
ray burst. 

It has also previously been suggested that the outflow can dis- 
sipate over a much wider range of radii, thereby creating a low- 
energy blackbody excess in gamma-ray burst spectra (Meszaros 
et al. 2002), or possibly a harder blackbody component just in- 
side the photosphere which could help to regulate the peak en- 
ergy in the observed manner (Meszaros & Rees 2000; Rees & 
Meszaros 2005). If the jet is magnetically dominated but very 
cold at its base, then continuing dissipation would be required 
to push it to a large Lorentz factor (Drenkhahn & Spruit 2002). 
However, it is difficult to establish clear predictions when the 
dissipation is driven by instabilities in the freely expanding por- 
tion of the outflow, given the number of physical effects which 
influence the growth of these instabilities. 

The energy density of a gamma-ray burst outflow is orders of 
magnitude higher than that of the Solar corona, and the ambi- 
ent radiation field has a much stronger influence on the cooling 
of charged particles. In fact, the conditions in the GRB out- 
flow are much closer to those expected in the magnetic corona 
of an accretion disk around a stellar-mass black hole. This sug- 
gests that, when drawing lessons from Solar physics, we should 
focus on those physical processes which appear to be energet- 
ically dominant in the Sun's atmosphere, and not only those 
which are directly responsible for the hard X-ray emission that 
emerges from that much more dilute environment. In the largest 
Solar flares, the output in thermal (<~ 1 keV) X-rays is smaller 
than the kinetic energy of the ejected plasma by an order of 
magnitude; and the output in hard X-rays and gamma-rays is 
smaller yet by six orders of magnitude (e.g. Somov 1992). Pro- 
cesses very similar to those which heat the Sun's chromosphere 
and thermal corona will, in the more extreme environment of a 
gamma-ray burst, create a much harder photon spectrum. 

A non-potential magnetic field transfers energy to particles 
through various channels. Energetically important channels in- 
clude the damping of large-scale magnetohydrodynamic mo- 
tions (through the formation of high-frequency spectrum of 
MHD waves); and a sudden change in the topology of magnetic 
field lines through reconnection. In a GRB outflow, charged 
particles are created at a sufficient rate by photon collisions 
(7 + 7 — > e + + e~) to ensure than the MHD condition applies on 
large scales. For this reason, we have previously suggested that 
the main role of reconnection in GRB outflows is not to cancel 
off the toroidal magnetic field (as in the Coroniti 1990 model of 
the Crab pulsar wind), but instead to facilitate a tangling up of 
the field through a change in topology. 

Damping of MHD turbulence is expected to be a ubiqui- 
tous process in any relativistic outflow containing a magnetic 
field. In this paper, we examine, in more detail than previously 
(Thompson & Blaes 1998; Lyutikov & Thompson 2005), the 
regime in which the magnetic energy density is at least compa- 
rable to the rest-energy density of the entrained charged parti- 
cles. The wave energy is transferred to the light charges (elec- 
trons and positrons) primarily through electrostatic acceleration 
along the magnetic field. As a result, the radiation is beamed 
along the local direction of the magnetic field. This means 
that inverse-Compton cooling will dominate synchrotron cool- 
ing even if the magnetic energy is larger than the energy den- 
sity of the ambient radiation field. Bulk relativistic motion of 
the magnetofluid provides another avenue for beamed emission 
(Lyutikov & Blandford 2003), but it is not required. 



An overarching goal of this paper is to identify the key com- 
ponents of the outflow which must be included for a proper 
description of GRB emission. We find that it is essential to 
include a thermal radiation field as well as a strong magnetic 
field, these two components being in approximate equipartition 
at the base of the outflow. A small admixture of baryons is 
probably present, which can carry a significant portion of the 
luminosity in some parts of the jet; but the baryons actually play 
a less essential role in the emission physics. For the purposes of 
investigating the prompt gamma-ray emission, the bulk of the 
outflow can be assumed to expand freely in between the photo- 
sphere of the Wolf-Rayet star and the dissipation zone. Closer 
to the central engine, the outflow probably is magnetically dom- 
inated and is susceptible to global MHD instabilities (Lyutikov 
& Blandford 2003; Giannios & Spruit 2006). The field lines are 
then entrained in the heated outflow, and stretched back into a 
non-radial configuration outside the Wolf-Rayet envelope. The 
fast motion of the jet implies strict limits on the factor by which 
the field energy can be reduced this close to the engine. 

1.1. Plan of the Paper 

The plan of this paper is as follows. We begin in §2 by 
defining the inner boundary conditions of the relativistic out- 
flow, and how they are influenced by the interaction between 
jet and Wolf-Rayet envelope. Basic constraints on the tangling 
of the jet magnetic field interior to the stellar photosphere are 
outlined in §2. 1 . The fundamental relation between black-body 
temperature an isotropic luminosity is worked out in §2.2. The 
baryons and the magnetic field are tightly coupled, and are ac- 
celerated together by the flux of thermal photons outside the 
baryonic photosphere (§2.3). Next we examine the interaction 
between the relativistic shell and the Wolf-Rayet wind (§3.1), 
with a particular focus on the radiative compactness in the de- 
celeration zone (§3.2) and optical depth in electron-positron be- 
hind the forward shock (§3.3). These calculations are repeated 
in §4 in the case where the fragmenting 'breakout' shell is the 
dominant source of inertia. 

The evolution of the magnetic field near the contact discon- 
tinuity is the subject of §5. We consider whether a reverse 
shock will form, and how its structure is modified by small- 
scale structure in the magnetic field (§5.1). We note that mag- 
netic reconnection across a neutral sheet can be facilitated by 
rapid synchrotron cooling in some circumstances (§5.3). The 
Rayleigh-Taylor instability of the contact discontinuity is an- 
alyzed in §5.4. The instability is weak except when a dense 
breakout shell is present. 

We argue in §6 that the damping of torsional MHD waves 
in a background radiation field provides a promising approach 
to the non-thermal and strongly variable emission of gamma- 
ray bursts. The section begins (§6.1) by reviewing the fact 
that electrons and positrons will be electrostatically accelerated 
along the background magnetic field by high-frequency waves 
near the inner scale of the turbulent spectrum (Thompson & 
Blaes 1998; Lyutikov & Thompson 2005). It is argued that 
electrostatic heating occurs primarily through many small im- 
pulses, at a scale somewhat larger than the scale at which the 
torsional waves become charge-starved. The beaming of the 
inverse-Compton radiation of the heated charges is calculated 
in §6.3. When the Kolmogorov energy flux is above a criti- 
cal value, the electrostatic heating cannot be balanced instanta- 
neously by Compton cooling (§6.4). 

The beaming of the inverse-Compton emission provides a 
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nice explanation for the rapid variability of the gamma-ray 
flux at higher energies. The correlation between variability 
timescale and photon energy is analyzed in §7.1. The effect 
of inverse-Compton scattering on pair creation in the relativis- 
tic ejecta is considered in §7.2, and it is noted that strong 
anisotropy in the gamma-ray emission can cause a significant 
increase in the pair creation rate. 

The paper concludes in §8 with a critical summary of how the 
model fares in explaining the non-thermal spectra of gamma- 
ray bursts. We show in §8. 1 that an inhomogeneous distribution 
of particle energies will naturally arise behind the reverse shock 
wave, and can explain the presence of a high-energy contin- 
uum above the peak energy in the spectrum of a typical GRB. 
The pairs swept up by the foward shock during the decelera- 
tion phase leave a radiative signature in the form of inverse- 
Compton photons peaking at an energy <~ 10 2 - 10 3 times the 
peak burst energy (§8.2). High-energy spectral measurements 
will therefore provide a direct diagnostic of the process of pair- 
loading and pre-acceleration. The implications of weak optical 
emission in the prompt GRB for the mechanism of particle heat- 
ing are drawn in §8.4, the transition from prompt GRB emission 
to afterglow is addressed in §8.5 and, finally, the 'peculiar' burst 
941017 is discussed in light of these results in §8.6. 

A summary list of our results is provided in §9. 

2. GEOMETRY AND INNER BOUNDARY CONDITIONS OF THE 
RELATIVISTIC OUTFLOW 

Our approach to the gamma-ray emission problem starts with 
a basic simplification: the relativistic outflow that powers the 
GRB suffers only a negligible dissipation between the central 
engine and the non-thermal emission zone at ~ 10 14 — 10 15 cm. 
The base of this freely expanding outflow coincides with the 
photosphere of the progenitor star. The temperature in the out- 
flow is already small enough that thermal pairs have annihi- 
lated (Goodman 1986; Shemi & Piran 1990), and the entrained 
baryons and electrons provide the main scattering opacity. We 
now discuss the relative partitioning of energy between the ther- 
mal photons, magnetic field, and baryons. 

In this paper we focus on a Wolf-Rayet progenitor, with a 
radius of <~ 2 x 10 10 cm. The star is a source of a strong stel- 
lar wind, M w ~ 0.5 x 10" 5 - IO^Mq yr" 1 (Nugis & Lamers 
2000) that provides an effective stopping agent for the relativis- 
tic ejecta. The central engine is a <~ 2-3M Q black hole that is 
surrounded by a neutronized torus (Woosley 1993). This torus 
is fed by the collapse of the inner core of the Wolf-Rayet star. 
Accretion and the emission of a Poynting-dominated jet contin- 
ues until the star explodes as the result of the energy deposited 
in it by the jet (Paczyhski 1998). The duration of the outflow is 
not lengthened compared with the period of activity Af of the 
central engine, and so we leave At as a free parameter, normal- 
ized to the observed durations of long GRBs, Af ~ 10 s in the 
cosmological rest frame. 

It should be noted that careful studies of the X- 
ray/optical/radio afterglow radiation of GRBs have called into 
question the presence of a dense Wolf-Rayet wind environment 
in several burst events. A lower mass-loss rate is, indeed, ex- 
pected from Wolf-Rayet stars of lower metallicity (Chevalier, 
Li, & Fransson 2004); but the mass-loss rates from luminous 
WC stars are also often inferred to be higher than ~ 1O~ 5 M 
yr" 1 . Observations of strong radio emission at <~ 1 month post- 
burst can place interesting constraints on the wind density (e.g. 
Chevalier et al. 2004), but under some assumptions, e.g. that 



the magnetic field energy fraction Sb is constant between the 
early and later synchrotron-emitting phases. (The wind density 
derived for GRB 020405 scales asM~ (sb,i month/£B,i day) 2 : 
see eqs.fl] and [2] of Chevalier et al. 2004.) 

The sharpness of the break in the optical light curve pro- 
vides a valuable probe of the density profile surrounding the 
burst source (Panaitescu & Kumar 2000, 2001, 2002). In the 
case of a jet with non-radial structure, it is difficult to repro- 
duce the sharpness of the break in about 1/2 of events, if the 
energy flux but not the Lorentz factor is allowed to vary as a 
power-law about the jet axis (Panaitescu 2005). On the other 
hand, the Lorentz factor r re i of the jet material is expected also 
to vary about the jet axis. Some calculations have included 
a strong angular gradient in r re i (Kumar & Granot 2003), but 
more general distributions of r re i have not been explored. Re- 
cently Swift has made available simultaneous X-ray and opti- 
cal light curves which have revealed different break patterns in 
these bands. GRB 050525a provides an example of an after- 
glow showing a much sharper break in the optical band than in 
the X-ray band ~ 10 4 s post-burst (Blustin et al. 2006). The 
implications for the density profile from jet break phenomenol- 
ogy therefore appear to us to be somewhat ambiguous (see also 
Yost et al. 2003). Given these caveats, we will focus solely 
on a r" 2 density profile in this paper and generally normalize 
M to a value ~ 10" 5 M Q yr" 1 . The implications of more dilute 
winds for the dissipation of the magnetized fireball are exam- 
ined briefly in §8.5. 

The simplest relativistic jet expands ballistically into a fixed 
solid angle. A linear increase of Lorentz factor with radius 
can be a good approximation to the actual situation even if the 
magnetic field is dynamically dominant close to the central en- 
gine. The transfer one half of the magnetic energy to thermal 
radiation allows the outflow to continue to accelerate by pres- 
sure gradient forces inside its photosphere (Drenkhahn & Spruit 
2002). We argue here that this process can be effective even 
when the thermalization is localized at some small radius, and 
that the acceleration can continue outside the outflow photo- 
sphere if the matter loading is very small. One then has 

r» = rw (i) 

up to some maximum Lorentz factor r max that depends on the 
thermalization radius. It is r max ~ 10 2 for a jet emerging from a 
Wolf-Rayet envelope. The remaining irregularities in the mag- 
netic field are then frozen in, and fall out of causal contact either 
because T saturates, or because the fluid passes through the re- 
verse shock wave (Thompson 1994; Lyubarsky & Kirk 2001). 

2.1. Tangling of a Dynamo-Generated Magnetic Field 

The most plausible mechanism for generating a baryon-poor 
outflow invokes a large-scale ordered magnetic field threading 
the black hole (Blandford & Znajek 1977), but such a Poynting- 
dominated outflow can undergo global current-driven instabil- 
ities as it propagates through the envelope of the Wolf-Rayet 
star. In our approach to the burst emission problem, these in- 
stabilities are not the direct source of the gamma-ray emission 
(as suggested by Lyutikov & Blandford 2003), but they will 
facilitate a tangling up of the magnetic field inside the Wolf- 
Rayet photosphere. The limitations of causality prevent the or- 
dered Poynting flux from decreasing to a small fraction of the 
total energy flux during the propagation of the jet through the 
Wolf-Rayet envelope. The field can be smoothed over a length- 
scale Ar' B (in the bulk frame) only if the Alfven crossing time 
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Ar' B jY' A is shorter than the flow time ?'. In a cold medium, the 
Alfven speed can be normalized to the speed of light through 
the fraction e' B of the energy density carried by the magnetic 
field: V'Jc ~ (2e' B )^ 2 /(l + 2e' B ) 1 ^ 2 . The duration of a flow that 
expands with Lorentz factor r is t' <~ r/cT, and so one obtains 
a lower bound e' B > (Ar' B /ct') 2 ~ (TAr' B /r) 2 on the magnetic 
energy density. If the field starts out with a reversal scale Ar' m 
and Ebq ~ 1, and if its reversals are random, then conservation 
of magnetic flux implies that e' B <~ Ar B0 /Ar' B . Combining these 
expressions gives 
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r/rj ■ (2) 

The magnetic field that is carried outward by the relativis- 
tic outflow is ultimately generated by a dynamo process within 
the central engine. The fluid shear in the neutron torus sur- 
rounding a 2- 3Mq black hole with a modest spin has a char- 
acteristic timescale of fd yn ~ 10" 3 s. The convection motions 
inside a nascent neutron star have a similar duration. The sign 
of the magnetic field that is advected into the outflow will vary 
with the sign of the dynamo generated field (Thompson 1994). 
The timescale for field reversals must be longer than the flow 
time r/c ~ 10~ 4 s at the base of the outflow; but simulations of 
the MRI instability also suggest that it will be much shorter 
than the total duration Af ~ 10 s of the inflow to the black 
hole from the collapsing core of the progenitor star (Woosley 
1993; Paczyhski 1998). This means that the characteristic ra- 
dial scale Ar B for toroidal field reversals in the outflow will not 
be much smaller than the radius of the Wolf-Rayet star. Tak- 
ing Ar B o/c ~0.1 s (before untangling) and R t = 2x 10 10 cm 
gives Ar B o/R+ ~ 0.15 and 
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from eq. (2). This lower bound is large enough that the field 
will have a significant effect on internal shocks farther out 
(§5.1). 

2.2. Seed Blackbody Photons and the Epeak ~ ^iso Af Relation 

A nearly black body photon gas is generated by the damp- 
ing of turbulent motions in the jet. This seed radiation field can 
carry a significant fraction e b b of the outflow luminosity L. At 
the stellar photosphere, r = R±, the rest-frame temperature T b ' b 
of the photons is given by 



T-i2 tI 4 

1 al hb = Ebb 



AnR 2 c 



: £bb 



4irR 2 c 



keV. 



(4) 



and the mean energy per photon carried by the thermal compo- 
nent of the outflow. This agrees closely with the relation be- 
tween £peak and isotropic burst energy L; so Af that was obtained 
by Amati et al. (2002) from a sample of BeppoS AX bursts with 
measured redshifts, and confirmed by HETE-II burst localiza- 
tions (Lamb et al. 2005). Second, the strong momentum flux 
of the black body photons allows this material to be accelerated 
to a higher Lorentz factor outside the Wolf-Rayet photosphere. 
The limiting Lorentz factor and the maximum tolerable matter 
loading are considered in the following section 2.3. 

Jet material moving at a Lorentz factor r> where 8 
is the angle from the axis of the jet, will sustain only a small 
velocity shear in the fluid rest frame, and will remain cold. Ma- 
terial that is decelerated to T <C 9~ l mixes easily with the much 
heavier cocoon material and decelerates to T <~ 1 . At an angle 
9, there is therefore a characteristic Lorentz factor T <~ 1 / \/36 
above which a Kelvin-Helmholtz mode will not have time to 
grow on a lengthscale ~ 6R+ (in the bulk frame). (Here we 
have equated the sound-crossing time 9R+/c s = V36R+/c with 
the flow time R+/Tc in the bulk frame.) 

If the entire jet fluid were to decelerate to T <~ 1 within 
the stellar envelope, then the relation between temperature and 
(isotropic) luminosity would follow the usual blackbody rela- 
tion, 7t,b ~ ^iso 4 - The temperature-luminosity relation becomes 
harder with T at the above critical value. Observations of GRB 
afterglows suggest that the total energy E carried by the jet is 
roughly independent of opening angle (Frail et al. 2001 ; Bloom 
et al. 2003). In that case, one has 

L iso Af (J^pj = L iso Af (^pj = 5 x 10 50 ergs. (5) 

The corresponding Lorentz factor is then 

r(^) = ^ = i-84 / o 2 5i^ 1 1/2 . 

Combining eqs. (4) and (5) gives 



(6) 



«=io4/ 4 V/ 4 



2x 10 10 cm 



-1/2 



keV. (7) 



This temperature is low enough that the thermal creation of 
electron-positron pairs can be neglected. The relation between 
Lorentz -boosted temperature and isotropic luminosity can then 
be expressed in terms of the mean energy per black-body pho- 
ton, 



2.7fc B r bb = 2.7 



= 70 



,1/4 
'bb 



(■ 



L iso Af 



Af 1/4 Vl0 52 ergs 



1/2 



2x 10 10 cm 



-1/2 



Here L; so = L (4tt / AO) is the isotropic luminosity of a two-sided 
jet flowing through a solid angle Afi = 2 x tt9 2 . Note that when 
the Lorentz boost of the jet material is taken into account, the 
temperature of the blackbody component could be as large as 
|r(/?*)7^ b ^100 keV. It therefore provides a useful seed for the 
prompt gamma-ray emission. 

As we now argue, the jet material that develops a high en- 
tropy inside the Wolf-Rayet star is favored as an observable 
source of rapidly variable GRB emission for two reasons. First, 
one obtains a simple relation between isotropic luminosity L; so 

1 The observed T90 distribution covers the range 4-40 s at half maximum (Paciesas 
source redshifts. 



keV. 



(8) 



If one substitutes Ebb ~ 5 and takes into account that the 
isotropic luminosities of long GRBs appear to have a much 
broader distribution than do their durations 1 , then one very 
nearly reproduces the observed relation (Amati et al. 2002; 
Lamb et al. 2005). The normalization is ^60% of the mea- 
sured value, but an additional <~ i of the outflow energy must 
be carried by a separate component, whose dissipation gener- 
ates the high-energy, non-thermal tail to the gamma-ray burst 
spectrum. We emphasize that this relation depends on only one 

et al. 1999), which of couse includes the broadening effects of a distribution of 
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free parameter, the radius of the Wolf-Rayet star, which is of 
course contrained independently. 

The inferred width of the burst energy distribution depends 
on the details of the beaming correction that is employed. The 
most detailed effort, by Ghirlanda et al. (2004), indicates a 
spread of <~ 1 order of magnitude in the case of a r~ 2 circum- 
burst density profile. However, Swift has observed additional 
time-structure in the afterglow light curves at times intermedi- 
ate between the burst duration Af ~ 10- 10 2 s, and the later 
achromatic breaks that were used by Ghirlanda et al. (2004) to 
normalize the jet opening angle (Tagliaferri et al. 2005). There 
remains, therefore, some ambiguity in the mapping between the 
prompt gamma-ray emission and the later phases of the after- 
glow. 

The distribution of burst energies also depends on the mech- 
anism by which the jet is collimated. The beaming-corrected 
burst energies are remarkably similar to the typical kinetic en- 
ergy of a core-collapse supernova (Frail et al. 2001). They are 
also similar to the gravitational binding energies of the cores of 
20-30 M Q evolved stars (typically £ s ~l-2x 10 51 ergs; see 
Figs. 5 and 6 of Woosley & Weaver 1995). Current modelling 
of Blandford-Znajek jets from black holes suggests that strong 
collimation requires the presence of a geometrically thick torus 
around the hole (e.g. McKinney & Gammie 2004). The extent 
to which this condition is satisifed in hyper- accreting tori (ac- 
cretion rates M ~ 2M Q /10 s = O.2M s" 1 ) depends on the den- 
sity profile in the outer atmosphere of the torus, which in turn 
is sensitive to the dissipative processes operating there. A geo- 
metrically thin, neutrino-cooled disk can form self-consistently 
at this accretion rate (e.g. Narayan, Piran, & Kumar 2001). 
When the outflow is broad close to the engine, strong colli- 
mation is possible only if the outflow energy is comparable to 
Eb\ if the outflow energy is smaller, it will be swept back by 
the accretion flow, and if it is substantially larger it will trigger 
a quasi-spherical blast wave. These considerations provide a 
physical motivation for our assumption of a relatively narrow 
distribution of burst energies, compared with current beaming 
corrections (Ghirlanda et al. 2004). 

Rees & Meszaros (2005) have outlined a related approach to 
the £peak ~~ ^iso Af relation, starting with the idea that the peak 
energy is regulated by dissipation and thermalization just inside 
the photosphere of the outflow. The radius of the photosphere is 
determined self-consistently by the luminosity and composition 
of the outflow, the scaling of T with radius, and the particular 
radiative processes occurring in the outflow. As a result, defi- 
nite predictions for the slope of the relation, or its normaliza- 
tion, are not easily obtained. Note also that the causal relation 
is more difficult to satisfy this far out in the out- 
flow. For all of these reasons, we conclude that the Amati et al. 
(2002) relation arises more naturally at a fixed radius, during 
jet breakout from the Wolf-Rayet star. The prompt gamma-ray 
emission is, in this picture, reprocessed 'fireball' radiation, with 
the position of the fireball displaced outward from the radius of 
the black hole accelerator (~ 10 7 cm) to its surrounding stellar 
envelope (<~ 10 10 cm). 

We expect that the relation (8) between T bh and L iso Af will 
be modified at low luminosities, because its derivation depends 
on r being larger than unity. One sees from eq. (6) that 
this assumption must break down below a luminosity L; so <~ 
At lower luminosities, the relation is 



An apparent inconsistency between the Amati et al. (2002) 
relation and the BATSE burst sample has been uncovered by 
Nakar & Piran (2005). These authors find that ~ \ of the 
BATSE bursts are harder spectrally than the £ pea k-£iso Af rela- 
tion would imply, no matter where each burst source is placed 
in redshift space. Note that the BATSE sample has a lower 
flux threshold than the BeppoSAX sample, so that the bulk 

1 12 

of the BATSE bursts need not sit in the £p ea k oc L iso portion 
of the color-luminosity relation. Our predicted modification 
of the Amati et al. (2002) relation, at energies lower than 
£j S0 ~ 3 x 10 51 ergs, could therefore resolve this apparent con- 
tradiction without recourse to beaming or orientation effects. 

The Tbb - Liso relation that we have derived is also consistent 
with a 'structured' jet, but we note that a range of Lorentz fac- 
tors within the jet core is required. The published structured-jet 
models generally assume a distribution of energy flux across the 
jet, but a constant Lorentz factor. In practice, there is no reason 
to expect that Y varies much more weakly than L[ so . The kine- 
matic model considered here implies that Y oc lV 2 A? 1/2 at the 
boundary of the Wolf-Rayet star. The asymptotic Lorentz fac- 
tor of the freely-expanding jet material scales as Y oc L 3 ^ 8 Af 
(see eq. [12] below). 

2.3. Baryon Contamination and Terminal Lorentz Factor 

Effective thermalization of the radiation field at a temper- 
ature (7) requires the presence of a small baryon component, 
since the temperature (7) is too small to excite pairs thermally. 
The mass in baryons depends, in turn, on the details of the 
launching of the jet, and can be expected to vary from burst to 
burst (depending especially on the viewing angle if additional 
material is mixed from the envelope). The baryon loading of 
the outflow is represented by the parameter 

M b c 2 

The associated electrons provide a scattering depth 

cr T Y e L iso ( r V 3 

TT ' h ^Y h Y\R*)m p R^ \Rj ' ^ 

Here Y e <~ | is the ratio of the number of protons to total nucle- 
ons that are bound up in charged ions. (The role of free neutrons 
is briefly addressed below.) 

A substantial optical depth at radius R+ is needed if the pho- 
ton distribution function is to relax to a black body. The jet 
material is still optically thin to free-free absorption at the tem- 
perature (7), and so effective thermalization requires a Comp- 
ton parameter yc = 4(kBT bb /m e c 2 )Tj y b > 5. Thermalization is 
therefore possible only if the baryon loading is larger than 



r 6 <7xio 4 4/ 4 Af: 5/4 



(11) 



3 x lO^Af" 1 ergs s" 1 



1/4 

predicted to soften to T^b ~ i iso when the jet material moves 
transrelativistically [Y(R+) ~ 1], 



at r = R± = 2 x 10 10 cm. The large scattering depth for ther- 
malization means that the blackbody photons must be created 
within the jet. More generally, unless Yb is extremely large 
and the jet is extremely clean, the jet core will not see the soft 
thermal photons that are created in the surrounding cocoon (cf. 
Ghisellini et al. 2000). 

Above a critical baryon density, the Lorentz -boosted photon 
temperature ^YT h ' b will be cooled by adiabatic expansion be- 
fore the photons and matter decouple. To avoid this, one re- 
quires that the scattering depth (10) drop below unity before Y 
reaches Y h . Combining eqs. (1), (6), and (10) gives the required 
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lower bound on I^, 



r h >r, 



frcrit : 



8irm p c 3 R+ 



1/4 



-1/4 



(12) 



The key point here is that Tbcnt depends on the launching radius 
of the freely expanding jet. Increasing R+ from ~ 10 6 - 10 7 cm 
(the physical size of the central engine) to R+ <~ 10 10 cm has the 
effect of reducing Tbcnt by an order of magnitude. 

The magnetic field carries inertia, which can limit the bulk 
motion of the outflow if it is cold at its base. The role of ther- 
mal pressure in accelerating the outflow has been examined by 
Drenkhahn & Spruit (2002). At a large optical depth, the radi- 
ation field is advected with the other components of the flow, 
and a pressure-gradient force is responsible for pushing the 
flow to a higher Lorentz factor. Although some acceleration 
is possible by this mechanism in the optically thin region of 
the flow, we focus here on the force imparted by the radiation 
field to the magnetofluid. This effect is already familiar from 
studies of particle-dominated fireballs, and is especially impor- 
tant if the outflow undergoes strong heating near its base. It 
is straightforward to see that the limiting Lorentz factor (12) is 
not modified significantly by the magnetic field as long as the 
particle Lorentz factor T is much larger than 2 (B 2 / '8tt pc 2 ) 1 7/3 . 
In a steady state, both the Poynting flux S = (B 2 /8ir)v r and 
the momentum flux P = (B 2 /8ir)(l + v 2 /c 2 ) carried by the non- 
radial magnetic field B receive (negative) corrections AS/S ~ 
AP/P ~ -Ar/r 3 in response to a change AT that is driven 
by Compton scattering. When T satisfies the above inequal- 
ity, both AS and AP are negligible compared with the increase 
pc 3 Ar in the particle kinetic energy flux. Further out in the 
flow, the inertia of electron-positron pairs can dominate the in- 
ertia of the baryons, and some of the energy carried by the mag- 
netic field is converted to non-thermal photons. 

It is possible to set an upper bound to the baryon kinetic lu- 
minosity, and to the terminal Lorentz factor of outflow. Outside 
the photosphere of the Wolf-Rayet star, the baryons will con- 
tinue to accelerate off the flux of thermal photons. The colli- 
mation of the photons grows with radius, T 7 (r) ~ T(R+) (r/R+). 
The baryons maintain a Lorentz factor T ~ T 7 as long as 



p < ]_ p CT (gblAso) 



2r 2 c 



(13) 



This leads to a maximum Lorentz factor essentially equal to eq. 
(12) (e.g. Shemi & Piran 1990). When T b > T bctit , the Lorentz 
factor saturates at a value 

r~r icrit (r b » i\ cri t), (14) 

and the baryons carry only a fraction e b ~ T b cnt/^b of the total 
energy of the outflow. The proportions of the energy flux car- 
ried by thermal photons, baryons, and magnetic field must of 
course sum to unity, 



Ebb + Sb + £B = 1 • 



(15) 



flux at the base of the jet, and the development of a large termi- 
nal Lorentz factor outside the Wolf-Rayet envelope. Portions of 
the jet which dissipate less inside the Wolf-Rayet star may at- 
tain a lower terminal Lorentz factor, and their deceleration will 
be pushed to a larger radius. 

Neutrons play a limited role in this jet model. Advected out- 
ward from the central engine, they are entrained by the protons 
and alpha particles at the Wolf-Rayet photosphere, and decou- 
ple from them somewhat inside the baryonic scattering pho- 
tosphere. The neutron mass flux M n can be much larger than 
the proton mass flux, due to the effects of neutronization in the 
central engine (e.g. Beloborodov 2003). We take a (represen- 
tative) ratio M n /Mp = 8, which corresponds to M n /M a = 4.5 
if all of the protons have combined into alpha particles. The 
neutrons then decouple from the alpha particles just after the 
combined neutron-ion fluid (temporarily) reaches the limiting 
Lorentz factor 



r„ = 



(M n +M a )c 2 



- 1 ^ T" 3 / 8 At 1 ' ' 



R* 



1/8 ^.1/4 

"bb i 2x 10 i cm 



-1/4 



T b 

T b cri 



In a cold flow, e B is related to the rest-frame magnetic energy 
density e' B by e B — 2e' B . 

An important feature of this acceleration mechanism is the 
correlation between the presence of a strong thermal photon 

2 Here p is the rest energy density of the entrained charges in the frame of the engine 



(16) 

The neutrons fall behind the outer half of the ejecta 
shell at a radius r 2 cAf = 5 x lO 13 /^ At\ ,A e^ 2 (7?*/2 x 
10 10 cm)" 1/2 (r A /r/ 7crit ) 2 cm. At this radius, a fraction 

ln2r„Af/f 1/2 = 0.l4 / o 8 51 At\ l% (R±/2 x 10 10 cm)" 1 / 4 of the neu- 
trons have decayed. The effective inertia of the decayed neu- 
trons (mass AM n ) is larger than the inertia of the alpha parti- 
cles, (T bcnt /T n )AM n > M a , when L iso > 10 51 erg/s. Even in 
this case, the neutron decay will have a modest effect on the 
dynamics of the contact discontinuity, which at this point has a 
significantly smaller Lorentz factor than that of the relativistic 
outflow (§4). The effect of the neutrons is also limited when Tb 
is larger than T bcn t and the inertia of the outflow is dominated 
by an advected magnetic field. 

3. DECELERATION OF THE RELATIVISTIC OUTFLOW 

We now consider the interaction between the relativistic 
ejecta and the Wolf-Rayet wind. The ejecta form a flattened 
shell outside a distance cAt from the star. A thin layer of 
shocked wind material collects in front of the ejecta shell, and 
is preceded by a forward shock wave. We focus first on the ide- 
alized, one-dimensional treatment of this problem, and defer 
consideration of Rayleigh-Taylor instabilities to §4 and §5.4. 

There is an extensive treatment of the dynamics of the re- 
verse shock wave in the literature, assuming ideal shock jump 
conditions (Sari & Piran 1995; Meszaros & Rees 1997; Sari 
& Piran 1999; Nakar & Piran 2004), and taking into account 
the effects of magnetic pressure in that approximation (Zhang 
& Kobayashi 2005). Our goal here is to include i) the effects 
of pair loading and preacceleration of the medium ahead of the 
forward shock; and ii) the possibility of strong dissipation of 
the non-radial magnetic field behind the reverse shock. The first 
effect is an inevitable consequence of the observed gamma-ray 
emission if it occurs inside ~ 10 16 cm; and the second is neces- 
sary if the emission is triggered by the interaction of the ejecta 
with the external medium (cf. §8.5). The relativistic outflow 
will also interact with the thin shell of material that is collected 
by the jet head following breakout. We defer a discussion of the 
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associated effects to §4, after we have analyzed the interaction 
with the expanding Wolf-Rayet wind material. 

We distinguish between the luminosity carried by the par- 
ticles and magnetic field; and the total isotropic luminosity L iso , 
which also has a contribution from the seed thermal photons. 
They are related by 



r re i S> r c 3> 1 in eq. (20), which gives the equilibrium solution 



Lrei = (1 -£bb)£iso- 

The mass profile in the Wolf-Rayet wind is 

Pw(r) = 



4irr 2 V w 



(17) 



(18) 



when the mass loss rate M w = M w s x 1O~ 5 M0 yr" 1 and wind 
speed V w = V w g x 10 8 cm s" 1 are constant over the emission re- 
gion. The mass swept up from the wind inside a radius r is 



M w {r)-. 



M w r 



(19) 



A relativistic jet emerging from the photosphere of a Wolf- 
Rayet star will also collect a thin shell of matter from the stellar 
envelope, which can exceed M w (r) for sufficiently small r (§4). 

3.1. Dynamics of the Contact Discontinuity 

Pair creation by gamma-rays streaming into the external 
medium has a strong influence on the propagation of the for- 
ward shock (Thompson & Madau 2000; Meszaros, Ramirez- 
Ruiz, & Rees 2001; Beloborodov 2002). What has not been 
examined is the feedback of the pre-acceleration of the pair- 
loaded Wolf-Rayet wind on the dissipative processes occurring 
in the relativistic outflow. We show here that the the bulk of the 
dissipation is delayed to a radius where the relativistic motion 
of the ambient medium has nearly, but not completely, decayed 
away. As a result, the compactness of the radiation emitted be- 
hind the forward shock is regulated to a characteristic value that 
depends weakly on L m \. 

The structure of the flow near the contact is displayed in Fig. 
1 . If the wind material were held static before passing through 
the forward shock, then T c would be obtained by moving to the 
rest frame of the contact and balancing the momentum flux on 
both sides, 3 



T 2 c p w c 2 ■■ 



:r t 2 (i-/3 rel /? c ) 2 



rel 



i 



4nr 2 J T 2 C {\+(3 C ) 2 



V 4-irr 2 J 



(20) 

We will separate our treatment of the deceleration of the con- 
tact from the acceleration of the relativistic outflow closer to 
the engine, and focus on the regime r re i = (1 - /3 2 el )~ 1//2 » T c . 
It should be kept in mind that, in the specific case of a break- 
out jet from a Wolf-Rayet star, the seed thermal photons are 
able to push the baryons and entrained magnetic field to a lim- 
iting Lorentz factor (12) that is close to the values of T c de- 
rived below. Continuing dissipation of the magnetic field com- 
bined with pair creation could push the particles to yet higher 
speeds (e.g. Drenkhahn & Spruit 2002); but the observed cor- 
relation between peak energy and isotropic gamma-ray burst 
energy (Amati et al. 2002; Lamb et al. 2005) suggests that 
this dissipation is not energetically important in between the 
photosphere of the Wolf-Rayet star and the gamma-ray pho- 
tosphere at <~ 10 14 cm (eq. [48]). We therefore approximate 



r - r - 

1 c — x eq — 



LielV w 

4M w c 3 



1/4 



=35L; e / 1 4 51 v,: / 8 4 M; 



1/4 
'-5 • 



(21) 



Reverse 'Shock' 



Breakout 
Shell* 



Forward Shock 




Contact 

FIG. 1. — Structure of the relativistic outflow between the forward and re- 
verse shock waves. The ambient medium is loaded with electron-positron pairs 
that are created by side-scattering of the gamma-ray photons. These pairs and 
the ions from the Wolf-Rayet wind collect in a shell behind the forward shock. 
In the idealized one-dimensional flow problem, this shell is separated from the 
shocked shell of relativistic eject by a contact discontinuity. The magnetic field 
that is advected across the reverse shock wave undergoes forced reconnection 
and tangles up. The reverse shock is, in fact, not a true shock but has a finite 
thickness (§5.1). The breakout shell is a thin layer of material from the Wolf- 
Rayet star that is advected out by the jet head, and becomes transparent at a 
radius ~ 1() 14 cm (§4). 

This result is modified when the effects of pair loading and 
pre-acceleration of the Wolf-Rayet wind material are taken 
into account. We are interested in the regime where the 
wind material reaches a Lorentz factor r am b 3> 1 (but still 
less than the Lorentz factor T c of the contact). The Wolf- 
Rayet wind material is then swept up by the foward shock 
at a radius close to where it was sitting at the moment of 
the explosion. Each shell of mass AMj passing through 
the forward shock now deposits a momentum r c r am b(/3 am b - 
P c )(M±AM h )c ~ -(r c /2r am b)(A^±AMi)c in the frame of the 
contact. The left-hand side of eq. (20) is therefore replaced by 
(r 2 /2r amb )M ±( o w c 2 , where 

M±(p,n e+ ) = —^ + l (22) 
P 

corrects for the inertia of the electron positron pairs (density 
n e +) that have been created in the ambient medium (Thompson 
& Madau 2000; Beloborodov 2002). The shell Lorentz factor 
becomes 

r e = (2r amb ) 1/4 r eq . (23) 

When pre-acceleration is efficient, T am b <~ T c /2, one has 



■M 



( L re ]V w 



1/3 



(24) 



\4M w c 3 

It turns out that A4± <~ 1 except very near the beginning of the 
deceleration phase, and so we define a reference Lorentz factor 

■tfrV3 T/V3»>-l/3 



1 eq - 1 eq ^rel51 m - 



-5 ' 



(25) 



3 This expression holds if the relativistic outflow and the ambient medium are both cold and particle-dominated in their respective rest frames. This condition will be 
approximately satisfied in the ambient medium even after it is loaded with electron-positron pairs and pre-accelerated (Beloborodov 2002). The momentum balance 
at the reverse shock also remains similar when the relativistic outflow is magnetically dominated, as long as the Poynting flux is converted efficiently to radiation. 
This is plausible if the magnetic field reverses sign on a small scale (e.g. Thompson 1994; Drenkhahn & Spruit 2002). 
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with the dependence on A4± removed. 

The baryons provide a useful set of markers for the flow: af- 
ter pre-acceleration their number density is 

Pw/mp 



n h = 



1-/3, 



amb 



— Zi amb 



(r amb > i). 



In the frame of the contact one has 

Pw 



r c ^- (r amb »i). 



m„ 



(26) 



(27) 



The rate at which baryons are swept up by the contact does not 
depend on the speed of the ambient medium, as long as the am- 
bient material moves only a modest distance from its starting 
point at the moment of the explosion. (Pre-acceleration does 
still alter the value of T c itself and the ram pressure of the on- 
coming flow.) 

The outflow that powers the gamma-ray emission in a long 
burst is concentrated 4 within a time interval At <~ 10 s and 
therefore carries a limited energy E le \ = L le \At. The energy 
which is swept up within a distance r from the central engine is 



AE re i(r) = ^ 



2T Z c 



= 2T l eq M w (r)c z , 



(28) 



with nearly equal proportions being dissipated behind the for- 
ward and reverse shocks. We are interested in the dynamics 
of the outflow close enough to the engine that electrons and 
positrons cool rapidly. Nonetheless, the ions carry the bulk of 
the energy flux across the forward shock in the latter part of the 
prompt deceleration phase (§3.3), and do not cool effectively 
by simple incoherent processes. About j of the outflow en- 
ergy is therefore converted to radiation, and a fraction e+ ~ ^ is 
transported by radiation across the forward shock, 

AE+ = e+AErei. (29) 
Above a critical value of the compactness 

a T AE+ ( e+ \ <7 T L re i 



4irr 2 m e c 2 \2r 2 / 4irrm e c 3 ' 
the exterior medium is loaded with electron positron pairs and 
accelerated to a high Lorentz factor comparable to r eq (Thomp- 
son & Madau 2000). This critical compactness is 4rit ~ 100 
when the gamma-ray spectrum has a high energy tail that ex- 
tends well above <~ 1 MeV (Beloborodov 2002). A Lorentz 
factor T am b ~ 10 2 is attained at £+ ~ 10 3 . We make use of the 
approximate scaling 

r amb - ( T— ) < HV <^ HV). (31) 



{h) (lo2 ^ lo3) ' 



which agrees at the <~ 20% level with the results obtained by 
Beloborodov (2002) for incident spectra with photon energy in- 
dices -1.5 to -2. 

One can obtain the equilibrium Lorentz factor (24) in a sim- 
ple way by comparing the energy AE + of the gamma-rays pass- 
ing across the forward shock, with the inertia of the ambient 
medium. This implies a maximum Lorentz factor 

AE+ 

r amb (max)= iA ,,,,, , (32) 



M±M w (r)c 1 ' 



or equivalently, 

r amb (max) 



£+ 



L T e\V w 



2s. 



+ 



=A . (33) 



eq ,2. vi: J r 2 r eq M w c 3 M± v , ei|/ 

4 We suggest in §8.5 that the end of the prompt gamma-ray burst coincides with 
star. A broader relativistic wind may, therefore, continue to flow from the central 
onto the central black hole. 



Substituting this into eq. (23), one sees that Y c ~ r eq is the 
largest Lorentz factor which can be maintained by a balance 
between momentum deposition on the two sides of the contact. 

The first phase of deceleration corresponds to the passage of 
the reverse shock through the shell. We now show that this oc- 
curs before pre-acceleration becomes ineffective. The prompt 
emission of gamma rays is therefore directly tied to the end of 
pair loading of the ambient medium. 

The beginning of deceleration can be identified with the ra- 
dius /?decei- at which the ambient medium can no longer be 
pushed to a Lorentz factor r amb = jf eq before hitting the for- 
ward shock. The corresponding compactness is (eq. [31]) 



^+(^decel-) 




(34) 



Expressing 1+ in terms of r eq using eq. (30) gives 

R decel -=lAxlO u e + Lllf 51 M 5 J* 5 V;l /6 cm. (35) 
Outside this radius, 

r c (r) = f eq 



^de 



-1/4 



(r > tfdecel-)- 



(36) 



The prompt phase of deceleration is complete when 
/ (2T 2 )~ l dr = r/3r 2 = cAt, which occurs at a radius 



'3T 2 cAt\ 



2/3 



-^deceit / Ji eq txA ' \ -2/3.1/3 ■>_] T7 A .2/3 



(37) 



^decel- V ^decel- 

or, equivalently, at 

/? decel+ = 2.5xl0 15 e! /3 L;/ 1 2 5 >- 1 _ / 5 6 y i ; 8 /6 Aff /3 cm. (38) 

Notice that /?decei+ depends more strongly on the relativistic out- 
flow parameters (L re i, Af) than it does on the pre-burst wind 
parameters. At this point, the Lorentz factor of the contact has 
dropped to 

r e (* decel+ ) = 53ei /6 L^ l M^Xi 12 W /6 , (39) 



and the Lorentz factor of the ambient material to 

T amb = 2.7 s 2 J 3 M 2 J" yi /3 Atf 3 . (40) 



Pre-acceleration has the smallest effect on deceleration for 
low values of M w and large burst durations Af . One can com- 
pare eq. (38) with the deceleration radius of a relativistic shell 
propagating in a static external medium. Setting cAt = r/2T 2 
gives 



eq 



decel : 



^relAfX 

M w c 



1/2 



= lxl0 14 Ll£ 1 At l MSvl{ 2 cm. 



w-5 w 8 



(41) 

The relativistic motion of the external medium has died away 
before the reverse shock completes its passage through the shell 
if 

A h Ml l _ 5 V-l > I2e + . (42) 

(The coefficient the right side of this equation is obtained by 
substituting £ + = £ CT n andr = 2r 2 q cAf in eq. [30].) 

the termination of the thermal photon pulse, due to the explosion of the Wolf-Rayet 
engine after the GRB ends, but at a diminishing rate due to the decreasing mass flux 
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It has been suggested that, close enough to the central en- 
gine, collimation of the emitted gamma-ray shell (outside the 
emission radius) would allow the ambient medium to be accel- 
erated to a Lorentz factor higher than T c (Beloborodov 2002). 
This requires the compactness 1+ of the photons escaping across 
the forward shock to exceed the minimum compactness that is 
needed to push the ambient medium to a Lorentz factor > T c . 
The further collimation of the photons outside the emitting ra- 
dius would then allow the ambient material to 'surf to a higher 
Lorentz factor, so that it detaches from the ejecta. This phe- 
nomenon can, indeed, be expected if the source is not sur- 
rounded by a dense wind and the pair-loaded medium has a 
small optical depth. But in the present context, the emission 
of gamma-rays cannot be decoupled from the pre-acceleration 
process: the optical depth of the pairs swept up from the 
Wolf-Rayet wind is larger than unity at r < /?decei-> an d the 
gamma-ray photons escaping across the forward shock are in 
fact trapped in the pair cloud that they create (§3.3). 

3.2. Compactness and Speed of Cooling 

At this point, we can check our assumption that heat elec- 
trons and positrons will cool rapidly after passing through the 
forward and reverse shocks. Transforming to the frame of the 
contact, and taking the radiation to be isotropic in this frame, 
the energy radiated is AE' Ki = (3 /AT C )AE K \. The radiation com- 
pactness is then 



3.3. Scattering Depth and Energy of the Swept-up Pairs 

The scattering optical depth between the reverse and forward 
shocks receives a significant contribution from electrons and 
positrons that are swept up from the ambient medium. This 
depth t± am b is regulated to a value close to unity, for the same 
reason that the radiative compactness is regulated. 

It is useful first to relate the deceleration radius to the radius 
R T =i at which the electrons in the pre-burst wind (Y e ~ i) pro- 
vide unit scattering depth, 

W^ = 10xl0 n^ i cm _ 

4irm p V w 

When this inertia of the wind becomes dominated by electron- 
positron pairs, the scattering photosphere increases signifi- 
cantly, to 



R 



R T= i = 3.7 x 10 14 M w - 5 V~l cm (48) 



during the first stages of deceleration. Comparing this rescaled 
photospheric radius with the deceleration radius, one has 



^decel- _ j ^£+T e q 



R 



= 0.15e + 



( £rel5lVi 



vr8 



w—5 



1/6 



(49) 



f(r) _ ^A£ r » _ 
A-Kr 2 m e c 2 



(43) 



One sees that the pairs swept up from the ambient wind provide 
7"±amb ~ 3-10 at r <~ /?decei-- The optical depth is even larger 
closer to the central engine. Pre-acceleration of the ambient 
medium is in fact suppressed inside a radius ~ /?decei— 
Further out in the flow, 

-l 



T±ambW = (-A4±-l) 



It is related to the compactness 
the forward shock by 

e = 



R 



T=l 



of the radiation field outside Substituting for fl d ecei- using eq. (47), one has 



4r e£+ 



(44) 



T±amb(r) = 3.3 



M±-l 



(- 



1/6 



(50) 



(51) 



when the emission is isotropic in the bulk frame. The critical 
value of £+ for pair loading and pre-acceleration is significantly 
greater than unity (Beloborodov 2002), and so one sees that £' 
itself is not much greater than unity in the frame of the flow. 
After substituting for AE' Ki and then for T c using eq. (36), we 
have 



£+ \Lrel5lVwgJ V^decel- 

The further evolution of the scattering depth depends on the 
evolution of the pair loading factor A4± — 1. One has (Be- 
loborodov 2002) 



M- 



1-0.1 



4n, 



(52) 



i\r) = 0.5 



4 



S+Leq 



^decel- 



-1/4 



The compactness is given by eq. (34) at the beginning of decel- 
eration, which implies M± - 1 = 0(1). Further out in the flow 
one has 



M- 



•1 ~ 



= 4.7- 



M 1/6 

m w-5 



ll LV 1 ' 6 V^decel 



£ +^rel51^8 



-1/4 



^decel- 



-1/2 



(53) 



(r > /?decel-)- 



(45) 



(46) 



At the final deceleration radius Rdecei+, we have 

AW) = 23 e-:' 6 L^\ M^l V-J 12 Atf 6 , 
with e+ — j. 

We see that the compactness is regulated to a characteristic 
value that depends weakly on the Wolf-Rayet wind luminos- 
ity and ambient density, and changes very slowly with radius. 
The compactness is large enough to ensure fast cooling of even 
mildly relativistic electrons and positrons. (The compactness is 
even higher in the portion the relativistic outflow that interacts 
with the fragmenting 'breakout' shell; §4.) This has impor- 
tant implications for the prompt gamma-ray spectrum, which 
we turn to in the following section. 



The scattering depth accumulated from the freshly swept-up 
pairs therefore decreases with radius as r± am b oc r" 3 / 2 outside 
r = Rdecei-- It drops below unity before the end of the prompt 
deceleration phase, and reaches the value 

t± amb (fldecei + ) = 0.1 L^ 5 ] Mf/_ 3 5 yj /3 At F 1 , (54) 
when the reverse shock has passed through the ejecta shell. 

The optical depth that a gamma-ray photon sees ahead of the 
foward shock (denoted by r± ex ) is significantly smaller than 
t± amb, because the ambient material advances only a small dis- 
tance beyond its initial radius ro before being intercepted by the 
forward shock. The wind material accelerates as it drifts back 
through the photon shell. The shell of radiation moving ahead 
of the shock at radius r > /?decei- h as , from eqs. (30) and (34), a 
compactness 

!/2 , .-1/2 

r 1 (55) 

^decel- / 
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and a width Ar = r/3T 2 . Having intercepted a radiation shell 
of compactness I < £+, the Lorentz factor of the wind material 
is 

2 , / „ \ 2 / \-l 



ramb(^0= ( - 



^decel- 



(56) 



Starting at a radius ro, the wind material advances a distance 
r-r 



Ar 






It- 




Tamb 


15 I 





amb 



1 / 



-3/2 



30 V^decel- ' 

The scattering depth seen by a gamma-ray photon moving 
ahead of the forward shock is, then, 

1/6 / \ -3 



■ r 



T±amb = 0.l£. 



M, 



w-5 



^rel5lV vv 8 



^decel- 



(58) 

One sees that the photons moving ahead of the ejecta shell 
are trapped inside a radius ~ /?decei— This brings us to the 
regime assumed by Thompson & Madau (2000), where the 
mean inertia per scattering charge behind the forward shock is 
~ (0.1- l)m e . Since the compactness of the radiation emitted 
at r > Rdecsi- is too small to push the ambient medium up to the 
speed of the contact, we conclude that there is never any de- 
tachment between the pair-loaded wind material and the ejecta 
shell. 

Finally, let us calculate the characteristic Lorentz factor of 
the pairs which are swept back across the forward shock. This 

is ^ 

2r amb (r eq ) (^decel-) ^ ^ 

During the last stages of the prompt deceleration, one has 
T± <~ 10 (as may be seen by substituting eq. [37] into eq. [59]). 
The cooling time of these particles is substantially shorter than 
the flow time, by a factor 1 /£'T±. 

Additional pairs are created when gamma-ray photons are 
upscattered above the pair-creation threshold at the forward 
shock, but the effect is only modest. For example, photons of a 
bulk-frame energy ~ m e c 2 have an optical depth ~ 1 to pair cre- 
ation at radius /?decei+ (§7.2). If a power-law electron/positron 
tail, dN/dj e oc 7 J 2 , is generated at the forward shock over the 
range of energies T± < j e < 7 cmax , then the increase in optical 
depth is Ar± am b/r± am b ~ 2/hi7 emax . Pairs will also be cre- 
ated in the region behind the contact, as a direct consequence 
of the gamma-ray emission. Some important subtleties which 
influence the pair-creation rate, associated with beaming and 
the inertia of the breakout shell, are addressed in §7.2. 

4. BARYON SHEATH FROM JET BREAKOUT 

We have, so far, treated the interaction of the relativistic out- 
flow with the Wolf-Rayet wind under the assumption that the 
outflow emerges with only a minimal baryon contamination. 
The jet head must expand sub-relativistically deep in the enve- 
lope of the star (Matzner 2003), but a thin outer shell of stellar 
material can be accelerated to relativistic speeds. The existing 
analytic treatment of jet breakout (Waxman & Meszaros 2003) 
is one-dimensional and does not address the sideways slippage 
of the the shocked stellar material at the head of the jet. We now 
give a simple estimate of the residual column S + of stellar ma- 
terial that survives sideways slippage away from the jet head, 



and compare this with the column of material that is swept up 
from the Wolf-Rayet wind (Fig. 2). 




j_ r~(3 e) 



Breakout Shell 



FIG. 2. — As the jet emerges from the star, the jet head accelerates to a 
Lorentz factor T > 1/ \/39. A residual shell of Wolf-Rayet envelope material 
becomes stuck at the jet head, and cannot flow to the side and out of the path of 
the relativistic outflow. This 'breakout shell' is optically thick close to the star, 
and is subject to Rayleigh-Taylor instabilities when it begins to cool at ~ 10 14 
cm. 

The forward shock wave accelerates close to the speed of 
light as it approaches the photosphere of the star (Tan, Matzner, 
& McKee 2001). Beyond this point, the sound speed of the 
shocked stellar material is c s ~ c/V3. The material in this 
shocked shell can flow a distance ~ c s t' transverse to the jet 
axis in the time t' ~ r/Tc that it takes the jet head to cross a 
radius r at Lorentz factor V. The material will therefore be lost 
from the jet head if T < 1/ v36, where 9 is the opening angle. 
To simplify matters, we refer to the isotropic jet luminosity Li S0 , 
and do not separate out the energy fluxes in thermal radiation, 
magnetic field, and baryons. We denote by AM* = AitR 2 ^* the 
equivalent spherical mass of stellar material riding at the head 
of the jet. This shell intercepts an energy L mo (f r /T 2 c) from the 
jet, where the coefficient / must be determined self-consistently 
from the dependence of T on radius. Neglecting radiative losses 
(the shell is initially optically thick) its Lorentz factor is 



= / 



AM*c 3 ' 



(60) 



This gives T oc r 1//3 and / = f • The proportion of the relativis- 
tic outflow which has been swept up by the shell grows slowly 
with radius; its thickness is 



Ar = r 

2r 2 



(61) 



Setting T = 1 /V3& at r = R+, and making use of relation (6), we 
obtain 

AM*c 2 _ 3^ 2 R+ 3 _ 3/2 -5/2 ( R* 

(62) 

The shell continues to accelerate as it intercepts further energy 
from the relativistic jet, so that 



1 / r 



V39 \R+ 



1/3 



1/3 



(63) 
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The reverse shock reaches the back of the ejecta (Ar = cAt ) 
at a radius 

-2 



fldecei(shell) = 7xl0 14 L 3 



At 



2x 10 1 



(64) 



Referring to expression (38) for /?decei+ (the deceleration radius 
in the absence of a breakout shell), one sees that the inertia of 
the shell will speed up the passage of the reverse shock through 
the relativistic wind. Indeed, the shell mass AM* remains larger 
than the mass (19) of the swept-up Wolf-Rayet wind out to the 
radius 



3/2 



Mj_ 5 V w& 



R* 



2x 10 10 cm 



cm. (65) 



It should be noted that the normalization of /?decei is very sen- 
sitive to the normalization of eq. (63), scaling as /?decei oc 
[IX/?*)] 6 . If, for example, one replaces c s = c/V3 with c when 
evaluating the width of the jet, then /?decei grows by a factor 
-30. 

It is clear from the above that Rayleigh-Taylor instabilities of 
the breakout shell can strongly modulate the gamma-ray emis- 
sion at larger radii. A homogeneous shell would become trans- 
parent to photons at the radius 



R, 



=1 (shell) = 8xl0 13 y c 1 / 2 L- 1 o / 5 4 1 Af- 3/4 



R+ 



2x 10 10 cm 



1/2 



cm 



where Y e <TTAM+/4irr 2 m p = 1. Here Y e 



(66) 

is number of elec- 



trons per baryon in the outer layers of the Wolf-Rayet star. Be- 
fore the breakout shell becomes optically thin, the particles and 
photons will thermalize at a low temperature. The effective 
temperature at the radius R T= i is given by 



■ [r(^ T=1 )] 2 ar e ' ff 4 = 



4nRi 



(67) 



Making use of eq. (63) to evaluate T(R T 
ature 



r e ff=-r(/? T=1 )r e ' ff =i.84/ 1 5 2 1 Af 1 1/2 



= i), one finds a temper- 



2x 10 1() cm 



-1/3 



Lorentz factor of the shell, we equate expressions (36) and (63). 
The corresponding radius sits between /?decei- (the beginning of 
deceleration assuming no breakout shell; eq. [35]) and ^decei+ 
(eq. [38]). One obtains 



^7-3/7 a .3/14 1/7 ,>-l/14 T/ l/14 

: 66L L 4 51 Af/ e< M w ! 5 Vj s 



R* 



2x 10 10 cm 



at the radius 5 



R r m 



1 x 10 1 



3/7 v 3/14 



r 3/14 
J iso51 



At 6/1 M 3/14 



2x 10 10 cm 



4/7 



-1/7 



(69) 



cm. 



(70) 



The compactness of the outflow in the bulk frame can be ob- 
tained by substituting these expressions for r max and Rr m into 
eq. (43), 



e\Ri 



)-9L 1/U 



At^e. 



6I1 MI'\V-I n 



R f 



2 x 10 10 cm 



-1/7 



(71) 

Note the very weak dependence of the compactness on the en- 
ergy of the relativistic outflow. The compactness at smaller 
radii is 

-3 / \ -1 / \ -2 

(72) 



Rr m 



Rr m 



keV 
(68) 

in the frame of the central engine. 

The pressure of the breakout shell is dominated by photons 
while it is optically thick. The shell therefore becomes much 
thinner as it expands, and it becomes subject to a corrugation 
instability, as we detail in §5.4. The thickness of a plume 
which emerges at radius r has a characteristic angular width 
86{r) < 1 /T(r). The (isotropic) mass of baryons that is col- 
lected by such a newly formed plume is lower than AM*, but 
only by a numerical factor em- Relativistic fluid that catches 
up with the shell at later times will continue to flow through 
this opening, and the radius Rd ece \ will expand by a factor e~jfi. 
This process can be expected to repeat itself on a smaller angu- 
lar scale after V has grown by one e-folding (that is, after the 
breakout shell has expanded a decade or so in radius beyond the 
transparency radius [66]). 

The last generation of Rayleigh-Taylor plumes may there- 
fore be clean enough that the final passage of the reverse shock 
through the relativistic shell follows our previous analysis. The 
peak Lorentz factor of the contact does not, however, attain the 
equilibrium value L eq (eq. [24]) in this case. To find the peak 

5 It should be recalled that our derivation assumes that the relativistic outflow moves at r re ] S> T 
derived previously. This approximation has the effect of reducing r max by 10-20%. 



5. DYNAMICS OF THE MAGNETIC FIELD NEAR THE 
CONTACT DISCONTINUITY 

There are several reasons to expect that the magnetic field 
will contribute a significant fraction of the pressure between 
the contact discontinuity and the reverse shock. First, most of 
the particle pressure behind the contact is supplied by electrons 
and positrons (§7.2) which lose energy rapidly to cooling. Sec- 
ond, simulations of relativistic jets show that magnetic and ki- 
netic pressures remain in approximate equipartition out a large 
distance from the engine (McKinney 2005a,b). Third, the mag- 
netic pressure increases by an order of magnitude across a re- 
verse shock that is only weakly magnetized (Kennel & Coroniti 
1984; Zhang & Kobayashi 2005). One caveat here is that strong 
pulsations in the outflow can force a reduction in a at a large 
radius: after the individual pulses spread and merge, the ratio 
of L P to M b c 2 drops by a factor of the pulse duty cycle. 

We found that the energy density in the relativistic outflow is 
approximately independent of its luminosity at the point where 
the reverse shock has completed its passage through the ejecta 
shell (when the inertia of the swept up material is dominated 
by the Wolf-Rayet wind). This result follows from the scaling 

1/2 

^decei+ oc LJ. in equation (38). Setting the Poynting flux carried 
by a toroidal magnetic field B equal to a fraction e B of the wind 
energy flux gives 

Z? 2 Lrt>\ 

e B -r^r- (73) 



47T 



There is a characteristic magnetic field in the outflow, 



B= 1 x 10 5 



b B m w-5 

d /3 ^ 6 Ar?/ 3 



G, 



(74) 



w8 "M 

in the frame of the central engine. This field is independent of 
L re i, and depends only weakly on M w and V w . (The field in the 

In fact, r max is close to the maximum jet Lorentz factor (12) 
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bulk frame is weaker by a factor <~ T" 1 .) The synchrotron en- 
ergy of a relativistic electron gyrating with a lorentz factor "f e 
in the bulk frame is 



^sync 



,0 3 7 2 — = 3 P-V 
7e m c c llOV 



4^ 



eV (75) 



(as measured by the observer). 



5.1. Reverse Shock Wave in an Outflow with a Stochastic 
Magnetic Field 

We now examine the significant changes to the structure of 
the reverse shock that can result from an advected magnetic 
field. The fast-mode speed in the fluid upstream of the reverse 
shock can be approximated by the cold fluid formula, 



B' 



■4 



1/2 



(76) 



c [4np'c 2 + (B' 2 )] 1 / 2 

(The fast mode is nearly isotropic in this regime.) Here p' = 
p'^+n'j.meC 2 is the total particle rest energy density in the bulk 
frame, p\ is the proper rest density of ions, n'± = n' e _ + n' e+ is the 
proper density of electrons and positrons, and 

(B'f 



a= 8^ = — B (7?) 

As previously, e' B is the fraction of the energy in the magnetic 
field in the bulk frame. Transforming to a relativistically mov- 
ing frame, 6 the ratio of of Poynting flux to (cold) kinetic energy 
flux is 2cr. 

It is possible for Vp to exceed the speed of the relativistic 
outflow with respect to the contact, so that a reverse shock 
wave does not form (Lyutikov & Blandford 2003). The ex- 
ternal medium will, instead, transmit momentum to the ouflow 
through a gradual compressive disturbance. This happens if 



rel 



< 



1 



(1+2(7) 



1/2 



(78) 



2T e - [l-(V F ) 2 /c 2 ]V 2 
We have seen that r rel (eq. [12]) and T c (eqs. [36], [63], [69]) 
may be comparable in the deceleration zone at <~ 10 14 - 10 15 
cm. This means that a does not have to be much greater than 
unity for the reverse shock to be suppressed in an outflow with 
an ordered toroidal magnetic field. 

It should be emphasized that the wave speed (76) applies to 
disturbances propagating in a uniform magnetofluid. The re- 
sponse of the fluid to compression is altered substantially if the 
field undergoes high-frequency reversals, over a distance much 
smaller than cAt . Stochasticity in the dynamo operating in the 
central engine will produce flips in the sign of the field that 
drives the outflow. Some further tangling of this field can occur 
as the jet pushes out through the envelope of the Wolf-Rayet 
star (§2). The weakened magnetic field is swept outward by 
the jet beyond the stellar photosphere. Starting in a disordered 
state, the field remains disordered in the bulk frame as long as 
r re i continues to increase linearly with radius. That is, the two 
non-radial components B'^ and B' g are comparable in magnitude 
and 



1 rel 



Bl 



(79) 



The field becomes predominantly non-radial in the bulk frame 
after the fluid Lorentz factor stops growing (e.g. after r re i satu- 
rates at the value [12]). 



The deceleration of the outflow involves a compression of the 
fluid on a scale <~ cAt , much larger than the field reversal scale. 
The reversing field is susceptible to reconnection and tangling 
upon compression (Thompson 1994). After tangling, the fluid 
becomes more nearly isotropic and its sound speed closer to 
c/\/3. This makes a crucial difference to the propagation of 
signals in the outflow. In this situation, the equation of state 
varies smoothly over a macroscopic distance within the mag- 
netofluid. A localized jump in density can therefore appear if 
the flow speed is higher than <~ c/ V5 upstream of the contact. 

In the GRB emission model discussed in this paper, the 
plasma is pair-rich and the mean energy per particle is a modest 
multiple of m e c 2 . A fluid approach is therefore appropriate to 
the dynamics of the electromagnetic field and the entrained par- 
ticles. This contrasts with a pulsar wind, where the rest energy 
of the particles comprises a much smaller fraction of the total 
wind energy. In that case, the mean Larmor radius of the parti- 
cles in a neutral sheet can approach the sheet thickness, thereby 
facilitating reconnection. Lyubarsky (2003) has argued that the 
reversing component of the magnetic field will be erased at the 
reverse shock in a pulsar wind, so that the shock jump condi- 
tion is equivalent to that of a flow containing only the mean 
magnetic field. 

We now consider the thickness of such a mock shock, 
in the case where the sign of the non-radial magnetic field 
varies stochastically on a scale Ar' B (as the result of a dy- 
namo process in the central engine; §2.1). The mean flux 
density averaged over a radial scale Ar' >• Ar' B is (B') ~ 
{2Lp/T 2 r 2 c) l l 2 {Ar' / Ar' B )~ x l 2 . The magnetization parameter 
associated with the smoothed magnetic field is 



/\2 



(g) 
%TT P 'C 2 



Ar' 
AV B 



(80) 



In this situation, the radially-averaged magnetic field decays 
much more slowly with smoothing length than it would in a pul- 
sar wind, where the field structure is nearly periodic (e.g. Coro- 
niti 1990). A reduction in a from a value <~ 1 down to <~ 0.03 
takes place over a lengthscale Ar <~ 30Ar B = 3(Ar B /0.1 s) 
light seconds. The degree of compression at the reverse shock 
will therefore depend in a non-trivial way on the dynamo 
timescale in the engine, relative to the duration of the burst. The 
case where the reverse shock is spread out over a lengthscale 
comparable to the thickness of the prompt ejecta shell deserves 
examination, but will not be addressed in this paper. 

This model of a shock in a flow with a reversing magnetic 
field implies several modifications of the mechanism of first- 
order shock acceleration that is commonly employed in GRB 
emission models. First, the heating of even mildly relativistic 
particles is de-localized. A distribution of particle energies will 
result from the finite rate of creation of electron-positron pairs 
per unit volume (§8.1). Second, this means that rapid Compton 
cooling can significantly suppress the upper bound on the parti- 
cle energy behind the shock, in comparison with a fast heating 
mechanism such as shock acceleration or resonant absorption 
of ion cyclotron waves by positrons. And, third, the heated 
light charges are not isotropic, being electrostatically acceler- 
ated along the background magnetic field (§6). 

Even when a standard magnetosonic shock can form, its 
compression is significantly weaker at large a than it would 



6 A number of authors use a to denote the ratio of Poynting and kinetic energy fluxes in the frame of the star (e.g. Kennel & Coroniti 1984; Drenkhahn & Spruit 
2002), but we are more concerned here with the properties of the magnetofluid in the bulk frame. 
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be in a cold fluid without magnetic field. Strong compression 
is essential for the formation of a hard particle spectrum by 
first-order Fermi acceleration (Blandford & Eichler 1987). The 
Lorentz factor behind a static shock is r ps ~ (2a-) 1 I 2 when the 
magnetic field runs parallel to the shock. The field is only com- 
pressed by a factor 1 + (4cr)" 1 (Kennel & Coroniti 1984). The 
compression ratio is essentially unaltered at small a - but only 
if a < 0.05, that is, if the energy carried by the magnetic field is 
less than ~ 10% of the kinetic energy of the entrained baryons. 
In particular, the value of a estimated in eq. (2) is large enough 
that the magnetic field will have a significant softening effect 
on the particle spectrum generated by shock acceleration. 

Calculations of the synchrotron emission from shock- 
accelerated pairs (Meszaros & Rees 1997; Kobayashi et al. 
2005, and references therein) will, therefore, require significant 
modification if the magnetic field carries 10 percent or more of 
the energy flux. In the rest of the paper, we will generally refer 
to the reverse shock without qualifiers, but it should be under- 
stood that it is likely to behave quite differently than a standard 
shock in an ideal magnetofluid. 

5.2. Shocked and Pair-Enriched Wolf-Rayet Wind 

The medium in between the contact discontinuity and the for- 
ward shock contains a relatively weak magnetic field compared 
with the relativistic outflow. The field that is swept up from the 
Wolf-Rayet wind is predominantly toroidal and carries a frac- 
tion 

B 2 

of the kinetic energy density in the wind. The progenitors 
of GRBs must differ in some key respect from the large ma- 
jority of Wolf-Rayet stars. Rapid rotation has been invoked 
(Woosley 1993). If the entire star rotated rapidly, then a rela- 
tively strong magnetic field would be advected into the wind. In 
the regime where the magnetic stresses themselves play a negli- 
gible role in driving the wind, one has Eb, w — B 2 R 2 V^/ 'M W V^ = 
10" 2 B\ A 1 „/2) 2 y r 2 t+7 M- 1 _ 5 V~l Here B* = B* 4 x 10 4 G is the 
surface magnetic field and V roU = V mt ^ x 10 7 cm s" 1 is the sur- 
face rotation speed of the Wolf-Rayet progenitor. (The chosen 
normalization of B+ corresponds to a net flux ~ 10 25 G-cm 2 , 
comparable to that threading an ordinary radio pulsar but a few 
orders of magnitude less than what is ultimately needed to drive 
a Poynting-dominatedjet with a luminosity of ~ 10 51 ergs s" 1 .) 

During the last stages of prompt deceleration, the external 
medium develops a relativistic motion (59) with respect to the 
ejecta shell. The magnetic field upstream of the forward shock 
can be obtained by noting that both the toroidal flux density 
and the baryon density are compressed by the same factor (eq. 
[27]), 



B' = 



Pw/r- 



-B w — Y C B W . 



(82) 



This field is further compressed behind the forward shock by a 
factor up to 3 (when T± ^> 1). 

Slow cooling of the shocked ions would prevent much addi- 
tional compression of the fluid behind the forward shock. We 
argue in §8.2 that resonant absorption of ion cyclotron waves 
by positrons (Hoshino et al. 1992) is not likely to prevent ef- 
fective thermalization of the protons behind the shock, due to 
the relatively small fraction of the kinetic energy carried by the 
positrons. 



The energy density behind the forward shock is 



K = r±(r cPw c 2 ): 



(83) 



V 2r am b / 

and so the fraction of the energy density carried by the imme- 
diate post-shock field is 

2 



£r = 



Q B' W ) 2 /87T 

UL 



- 18r am b£B,w -2 X 10 4 r am b£B ! wV' 2 8 . 

(84) 

We have found r am b ~ 2 during the final stages of the prompt 
deceleration. The linearly compressed magnetic field can there- 
fore attain a pressure in excess of 10~ 5 of the thermal pressure 
behind the forward shock. This is close to some estimates of 
the magnetic energy density from the observed optical of some 
GRBs (e.g. Panaitescu & Kumar 2004, and references therein). 

The structure of the magnetic field in the Wolf-Rayet wind 
can be contrasted with that expected in the relativistic outflow. 
The Wolf-Rayet wind takes a year or so to flow out to the radius 
^decei+ ~ 10 15 cm (eq. [38]). Even if the progenitor supports an 
active dynamo, it will be located deep in the star. The magnetic 
field anchored at the base of the Wolf-Rayet wind is therefore 
likely to be fairly constant over the time the wind expands to 
the deceleration zone of the relativistic GRB ejecta. 

It has been noted that the Weibel instability will be excited 
behind a relativistic shock wave that is moving into a weakly 
magnetized medium (Kazimura et al. 1998; Medvedev & Loeb 
1999). This electromagnetic instability is generated by the 
two counterstreaming populations of electrons (and positrons) 
behind the shock, and creates strong current fluctuations on 
the plasma scale ~ c/ui Pe . The instability can be modified or 
suppressed by a strong seed magnetic field B' w = T C B W if the 
associated electron gyrofrequency uj ce = eB' w /m e cj e is larger 
than the plasma frequency ui Pe = (4im'±e 2 /m e ^ e ) 1 ' 2 (Hededal 
& Nishikawa 2005). Taking % = T± = r e /2r amb one finds 



LUp e 



B',„ 



47r(r c /2r amb )n' ± m e c 2 



1/2 



1/2 



,w 

M±-l 

(85) 

0.2 (the pairs carry 
<~ 2, this expression 



(86) 



Substituting (M± - 1) ~ n' ± m e /(T c p w ) ~ 
20 percent of the particle inertia) and r amb 
becomes 

^-0.02*^. 

LO Pe 

The pre-existing magnetic field is therefore not strong enough 
to suppress the Weibel instability. 

Nonetheless, the small grow scale ~ c/top e of the Weibel in- 
stability suggests that both the mean value and the r.m.s. value 
of the magnetic field will relax to very small values downstream 
of the shock, leaving the flux that is swept up from the Wolf- 
Rayet wind. Particle-in-cell simulations of the aftermath of the 
Weibel instability (Medvedev et al. 2005) have demonstrated 
that the coherence length of the magnetic field will grow with 
time; but these simulations cover only a limited dynamic range. 
On very long wavelengths, the smoothing out of the magnetic 
field is limited by the speed Va = B' / y/4-np' with which the field 
can move the entrained charged particles (Gruzinov 2001). To 
estimate the limiting magnetic field strength, we suppose that 
the field has reached equipartition with the kinetic energy den- 
sity of the pairs on the plasma scale 

/ 2 \ V2 
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namely 



B'(£ p ) ~ B' eq = {%n^ e ri±m e c 2 ) 



1/2 



(88) 



The conservation of magnetic flux implies that the smoothed 
field B' is 

-i 2 



B\L) 
B' eq 



<7jn±L \m e cL 



n 



(89) 



on a scale L > £ p . The associated Alfven speed is Va(L) = 
B\L)/(Airp') x l 2 . Setting L/V A (L) equal to the post-shock flow 
time <~ r/Y c c (as measured in the bulk frame) gives 

2 

T±m e c 



B'(L) 



eq 



(M ± -iy/ 2 eB' eq (r/T c ) 
(r c /50) 3 



= 5 x 10 



r l5 {M±-l)V 2 \ 2 



Vl0 5 G 



(90) 

This shows that the smoothed field that is left behind from the 
Weibel instability is weaker than the field that is swept up from 
the Wolf-Rayet wind (eq. [84]), even if the rate of smoothing is 
only <~ 10" 4 of the maximal rate V A /L. 

5.3. Feedback of Particle Cooling on Magnetic Reconnection 

The type of dissipation that is triggered by magnetic recon- 
nection requires careful clarification. In the equatorial region 
of a pulsar wind, it has been suggested that the reversing mag- 
netic field will largely cancel out (Coroniti 1990). In this man- 
ner, the fluctuating component of the non-radial field may be 
erased whilst preserving the mean field (Rees & Gunn 1974; 
Lyubarksy 2003). The alternative, which we favor in a rela- 
tivistic outflow with a much larger particle density, is that re- 
connection changes the topology of the magnetic field lines, so 
that the non-radial field is partly converted to radial field. In an 
expanding outflow, this allows the magnetic energy to decrease 
more rapidly than it otherwise would (Thompson 1994). 

A minimal particle density n em i n = \J r \/ec is required to sup- 
ply the radial electric current in an outflow carrying a toroidal 
magnetic field, 

(A:elC) 1/2 



\Jr\ = ' 



Airr 



dB d 
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(91) 



The actual pair density in a pulsar wind may exceed n em i n by 
4-5 orders of magnitude (e.g. Hibschman & Arons 2001). In 
a gamma-ray burst outflow there is generally a radius at which 
the Thomson depth is unity, outside of which the electrons (and 
pairs) are largely frozen into the flow. At this photosphere (ra- 
dius R T= i) one has n e = (cr T ,r? T= i/2r 2 )~ 1 , and 

n e %Tte[T(R T=l )] 2 R T=l 



= 2x 10' 



(Tt(Wc) 1 / 2 

6 [W=l)] 2 



r l/2 



2x 10 10 cm 



(92) 

In a gamma-ray burst outflow, one can also expect B^ to 
switch sign on a scale Ar B <C r, due to the intrinsic stochasticity 
of the dynamo in the central engine (§5.1). The current then has 
a dominant component \ Jg\' <~ cB'^/ At: Ar' B = cB / 'AnT 2 Ar B in 
the bulk frame. Comparing \Jg\ with \J r \ = \J r \/T, one sees 
that the minimum charge density needed to support the current 



increases by a factor <~ r/TAr B . For realistic values of Ar B , 
this factor is ~ 3 x lOV^T^ 1 [(Ar B /c)/Q.\ s]" 1 , which is sig- 
nificantly smaller than the pre-factor in eq. (92). The charges 
advected by the outflow are therefore capable of maintaining 
the non-radial current out to a radius where photon collisions 
will raise the pair density above the freezeout value (§7.2). 

The importance of pair creation has sometimes been for- 
gotten when considering macroscopic breakdowns of MHD in 
gamma-ray burst outflows. For example, strong wave accel- 
eration at the forward shock (Smolsky & Usov 1996) will be 
suppressed by pair creation outside the forward shock; and sim- 
ilarly long-wavelength electromagnetic disturbances (Lyutikov 
& Blackman 2001) will behave in an essentially hydromagnetic 
manner. The electromagnetic field cannot develop an essen- 
tially electric character (E 2 > B 2 ) because the entrained electric 
charges would have to move at a velocity exceeding the speed 
of light. Nonetheless, a localized breakdown of the MHD ap- 
proximation is possible where the magnetic field in the outflow 
has discontinuities. 

We now consider the effect that radiative cooling has on the 
speed with which magnetic flux can be advected toward a neu- 
tral sheet (given that such a field configuration is present in a 
GRB outflow). A pure electron-positron plasma trapped in a 
neutral sheet can cool rapidly, which has the effect of reducing 
the volume integral of the particle pressure across the sheet and 
greatly increasing the compressibility of the medium. By con- 
trast, ions will cool very slowly by incoherent synchrotron or 
inverse-Compton processes in the emission region of a gamma- 
ray burst. The ions can therefore supply a significant back pres- 
sure (Fig. 3). 

The relativistic expansion of the outflow has an important 
influence on the geometry of the reconnecting magnetic field. 
In contrast with the Solar magnetosphere, the magnetic field 
lines are nearly homogeneous away from the neutral sheet, and 
extend far beyond the causal scale r/T c . As a result, one ex- 
pects multiple X-type points to form along the sheet (Thomp- 
son 1994; Lyubarsky 2003). The rate of particle cooling implies 
a characteristic thickness A for the layer within which the mag- 
netic field becomes disordered. 




FIG. 3. — Distribution of magnetic field lines around an idealized neutral 
sheet. The magnetic field is turbulent within the sheet, but all three compo- 
nents of the spatially-averaged field vanish at the mid-plane. When the en- 
trained charges cool rapidly, it is possible for the oppositely directed field lines 
to reconnect and cancel out, without any sideways fluid motion. The thick- 
ness A of the sheet is determined self-consistently by the balance between 
the Poynting flux toward the sheet, and the outward radiative flux. However, 
a small contamination by slowly cooling ions will induce a significant back 
pressure which chokes off the transport of magnetic flux toward the neutral 
sheet. 
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In what follows, we assume that the magnetic field becomes 
strongly turbulent within this layer, and that the energy in the 
fluctuating component of the field is transfered to the particles 
by a Kolmogorov-type cascade on a timescale <~ A /c. As a re- 
sult, the magnetic field has annihilated in a layer of thickness 
<~ ct ' after the lapse of a time t ' < r/cT. The buildup of particles 
in the sheet forces a decrease in the magnetization parameter 7 , 

(93) 

Here a± is the ambient value, as determined self-consistently 
by the gamma-ray spectrum. Balancing the advection time A/c 
with the synchrotron cooling time <~ (r/cY)(£' B a±)~ 1 , where 



gives 



and 



e _ cn(B') 2 r 
B 8-Km e c 2 r, 

A 2 1 



<7± 



(r/T)(ct>) £' B a° ± 

(4) 1/2 (sL 

(4)V2 [r/r 



-1/2 



(94) 
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(96) 



The scattering depth across the sheet is smaller than the total 
depth of the outflow (t t ~ £' B /a±) by a factor Yet' / r, 



r T (A) . 



4 



ct' 

7Jf 



(97) 



The optical depth to annihilation (over a timescale t ') is 

01 \2 / ./ \ 2 



n(f) = 



(a ann v) r T (A) (£' B ) 2 



4 



ct 

rjY 



■ (98) 



a T c A/ct' 

Here we have made use of (<r ann v) <~ ajc/a± for relativistic 
pairs of energy <~ a±m e c 2 . One sees that the annihilation of the 
collected pairs can be neglected as long as a± 3> (£' B ) 2 - 

Now consider the effects of a small contamination by ions, 
with a magnetization parameter a p = B' 2 /4-nn' p m p c 2 3> 1. The 
ions collect in the neutral sheet. Reconnection can be choked 
off if the ions acquire a significant fraction of the wave pres- 
sure. Cooling of the ions by an incoherent process like syn- 
chrotron emission is very slow, and can be neglected as long 
as Up < (m p /m e ) 3 £' B . The dominant heating mechanism of the 
pairs is likely to be electrostatic when cr±> I (§6). The ions 
will be accelerated to relativistic energies by the same electric 
field if a± > m p /m e in the neutral sheet, but will absorb a sig- 
nificant fraction of the wave energy only if the numbers of ions 
and positrons are comparable, i.e., only if a p > (m p /m e )4- 
It should be noted that requiring a± > m p /m e translates into 
a stringent constraint on the magnetization of the pairs out- 
side the sheet, 4 ^ (m p /m e ) 2 £' B (from eq. [96]). More ef- 
fective ion heating would result from resonant absorption of 
high-frequency fast waves (e.g. Yan & Lazarian 2002). 

Finally, it should be noted that two-dimensional models of 
reconnection at a neutral sheet oversimplify the field geometry 
that is expected in a GRB outflow. If the magnetic field starts 
in a disordered state, then the field at the boundary of a fluid el- 
ement in the jet will typically be inclined by a large angle with 
respect to the field in a neighboring fluid element. (In other 
words, the two non-radial components B\ and B' g are compara- 
ble in magnitude, in a spherical coordinate system aligned with 

7 Throughout this paper, we define the magnetization parameter in terms of the 



the jet axis.) Near the boundary between two fluid elements, 
this means that one tangential component of the field will keep 
a constant sign and will not vanish in the current sheet at the 
boundary. The uniform advection of magnetic flux toward the 
fluid boundary will therefore be suppressed by the pressure of 
this component of the field. 

5.4. Rayleigh-Taylor Instability 

We now consider the Rayleigh-Taylor stability of the contact 
discontinuity that sits at the front of the relativistic ejecta shell. 
We consider two separate cases: where the inertia of the mate- 
rial that is swept up in front of the contact is dominated by the 
Wolf-Rayet wind (§3.1), and where it is dominated by a 'break- 
out shell' that is collected during the emergence of the jet from 
the outer layers of the star (§4). We then consider the influence 
of a magnetic field on the instability. 

The effective gravity felt by material near the contact discon- 
tinuity is 



d 2 r' 
(dt') 2 



> dY c 
dr 



(99) 



where dr' is the radial displacement in the bulk frame. This 
is directed outward when the effect of the breakout shell is 
negligible and the dynamics of the contact is dominated by 
the interaction with the Wolf-Rayet wind. As radiative pre- 
acceleration is turning off, the Lorentz factor of the contact 
decreases with radius, as Y c oc r~ x l A . The rest-frame time co- 
ordinate is t' = J (dr/Y c c) ~ |r/T c c, and so 



5t' 



(100) 



On the other hand, the Lorentz factor is still growing close to 
the central engine, r e oc r 1/3 , where the breakout shell domi- 
nates the inertia of the contact. The effective gravity is inward 
in this regime, and also somewhat stronger, 



c 

~2t' 



(101) 



Here we have made use of t ' ■ 



r/T c c. 



The forward shock is moderately relativistic near the end of 
the prompt deceleration phase (eq. [59]). The pressure P{ on 
the front side of the contact discontinuity is dominated by slow- 
cooling ions of an energy Y± ~ 10. The gravitating mass den- 
sity p' + can be expressed in terms of P' + and the energy density 
e' + ~ 3Pi by 

p ' + ^L 2 (e' + + Pi)^^±. (102) 

On the back side of the contact, the particles (mainly pairs) are 
only mildly relativistic when the effects of Compton cooling are 
taken into account. Their pressure is P'_ < \p'_c 2 in that case. 
The two sides of the contact will be in approximate pressure 
equilibrium, P'_ = P' + , and so one deduces that the fluid behind 
the contact is denser, p'_> p' + . The effective gravity is outward 
during the last stages of pair-loading in the external medium, 
and so one expects the contact to be Rayleigh-Taylor unstable. 

The fluid on either side of the contact suffers from a strong 
corrugation instability when its thickness Ar' is much smaller 
than r/Y c in the bulk frame. This instability is demonstrated 
by Vishniac (1983) in the case of a thin cooling shell in a non- 
relativistic blast wave, but the derivation is not altered in the 
relativistic case when Ar' <C r/Y c . In this regime, the growth 

rest energy density of the particles in the bulk frame, a± = B 2 /4irn ' ± m e c 2 . 
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rate is determined by the net acceleration of the shell material 
(see eqs. [2.19]-[2.21] of Vishniac 1983). 

The peak growth rate is 7' <~ {g'k\\) x l 2 ~ (2irg' / Ar') l l 2 at a 
wavenumber ku ~ 2ir/Ar' parallel to the contact. When the 
contact is decelerating, the ratio of flow time to growth time 
is 7V ~ (27r/5) 1//2 (Arycf')" 1/2 - This shows that the growing 
Rayleigh-Taylor fingers will induce only a limited amount of 
mixing when both fluid shells (on either side of the contact) 
have a thickness Ar' <~ ct'. The corrugation instability of an 
adiabatic relativistic shell has been analyzed by Wang, Loeb, 
and Waxman (2002), who find no growing modes. But growing 
modes do definitely exist in the thin-shell limit. 

The breakout shell, in particular, suffers from a strong cor- 
rugation instability. The direction of gravity and the density 
contrast are g' < 0, p' + > p'_ in this case, since the shell material 
is cooled by adiabatic expansion. The shell also becomes rela- 
tively thin as it is pushed outward. The pressure inside the shell 
decreases with radius as P' ~ L K \/AirT 2 c r 2 c oc r -8 / 3 , given the 
scaling Y c oc r 1 / 3 (eq. 63). Before the shell becomes optically 
thin, P' is dominated by the entrained photons. The pressure 
therefore scales with the shell thickness Ar' as P' oc ( Ar'r 2 ) -4 / 3 
under adiabatic expansion, and the aspect ratio decreases as 
Ar' /ct' oc r" 2 ' 3 . The shell suffers from a strong corrugation 
instability even before it becomes optically thin to scattering 
(at radius [66]). 

The details of the Rayleigh-Taylor instability change when 
the magnetic field dominates the pressure behind the contact. 
If the field were predominantly toroidal and of a constant sign 
(running parallel to the contact), then a Rayleigh-Taylor insta- 
bility would be suppressed. This can be seen from the disper- 
sion relation (Chandrasekhar 1981) 



is then 



L0' 2 = g'k ll 



P+-P- 



+2v 2 _q 



(103) 



K p' + + P'-J " 11 \p' + + pL. 

which is the sum of the usual Rayleigh-Taylor piece and a 
second positive piece representing the restoring force of the 
bent magnetic field lines behind the contact. The growth rate 
7' = Imci/ takes the maximum value 



Tmax ' 



at a wavevector 



P\~P- 
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Effective growth of the Rayleigh-Taylor mode requires 
7'f ' 3> 1 . When the contact is decelerating (g' is given by eq. 
[100]), this is possible if V A -/c < 0.1 \{p' + - p'_) / p'_\. The par- 
ticles behind the contact are mildly relativistic in this situation, 
and so a large density ratio p'_ 3> p' + is inconsistent with pres- 
sure equilibrium across the contact. We conclude that a modest 
magnetic field, which supplies ~ 10~ 2 of the pressure behind 
the contact, can suppress the Rayleigh-Taylor instability. 

On the other hand, a reversing magnetic field is susceptible 
to reconnection and tangling after it passes through the reverse 
shock wave. After the field has broken up into small-scale 
loops, it then behaves more like an isotropic fluid on larger 
scales. The stress-tensor of the magnetofluid must still be some- 
what anisotropic, with a non-radial component than is larger 
than 2 1 / 2 times the radial component, B'^ + B' g 2 = 2(l + e)B' r 2 . 

2 2 2 

We can define an average pressure P'_ = (B^ +B' e -B r 
normal to the contact. The equation of state of the magnetofluid 



B'/ B', 2 3 + 2g 

8tt 8tt 8tt l+2e 



(eL+PL). (106) 



Equating the pressures on either side of the contact and relating 
them to the densities through equations (102), (106), gives 

4-^. (107. 
PL l+e 

This expression assumes that e is larger than I/7 of the rela- 
tivistic particles forward of the contact, so that p' + ~ AP'+fc 2 . In 
this case, the anisotropic magnetic field is effectively the lighter 
fluid. The Rayleigh-Taylor mode is therefore stabilized while 
the contact is decelerating. 

6. DAMPING OF ALFVEN TURBULENCE IN A 
MAGNETICALLY DOMINATED MEDIUM 

Long-wavelength magnetohydrodynamic turbulence is a 
generic outcome of various instabilities in a GRB outflow, 
which include the tangling of the stretched non-radial field lines 
by reconnection (Thompson 1994) and kink instabilities in a 
jet (Lyutikov & Blandford 2003). The outflow retains enough 
charges outside its photosphere to enforce the MHD condition 
on very small scales compared with the causal scale (§5.3). We 
now examine how energy is transferred from waves to particles, 
and thence to the photons, when the magnetic energy density 
approaches or exceeds the rest energy density of the entrained 
charges. 

We should emphasize that we are examining the regime 
where the energy density in the background magnetic field 
is much larger than the energy density in high-frequency tor- 
sional waves. The conservation of Kolmogorov energy flux, 
from some forcing scale (eq. [126]) to higher frequencies, im- 
plies that the r.m.s. wave amplitude decrease with frequency 
as {5B 2 )H 2 /B Q oc lu 1 ' 2 . Models of 'jitter' radiation (Medvedev 
2000) focus on the opposite limit where the magnetic field has 
only weak long-range order, and the gyroradius of a particle is 
much larger than the coherence length of the field. 

Our focus here is on the torsional MHD waves which are ex- 
cited when a reversing magnetic field in the outflow is forced 
to reconnect. These waves transport energy along the back- 
ground magnetic field. The excited modes have a characteris- 
tic wavevector k±o <~ 1 / Ar' B perpendicular to the background 
magnetic field. Here Ar^ is the reversal scale of the magnetic 
field (due to the stochasticity of the dynamo operating in the 
central engine; §2). The component of the wavevector parallel 
to the background field is fcyo ~ (V le c/c)k±o, where V rec is the 
speed of transverse motions excited by reconnection. The ex- 
cited turbulence will therefore be anisotropic even at the outer 
scale Lq <~ l/ku . In what follows, this forcing scale will be 
normalized to the causal scale of the outflow, 



Lo = £0 =r ■ 

1- r 



(108) 



The amplitude of the turbulence at the forcing scale is SBq. No- 
tice that 



£0' 



Ar 



= 0.3- 



ri5 



Vrec/C 
0.1 



Aryc 
0.1 s 



(109) 

is characteristically ^0.1 if the magnetic field reverses on a 
timescale of ~ 0.1 s. In this section, we work consistently in 
the bulk frame, and so we drop the ' on all quantities for nota- 
tional convenience. 
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Torsional MHD waves can be damped directly by bulk 
Compton drag (Thompson 1994); or indirectly when non-linear 
couplings between oppositely propagating waves create higher 
wavenumber modes (Goldreich & Sridhar 1995). The current 
density increases with wavenumber in such a cascade, and an 
upper bound to the frequency of the waves is obtained by bal- 
ancing this fluctuating current with the maximum conduction 
current available from the embedded charges (Thompson & 
Blaes 1998; Lyutikov & Thompson 2005). Damping of the 
wave motions occurs primarily through electrostatic accelera- 
tion of these charges along the background magneic field. The 
accelerated charges cool primarily by inverse-Compton scatter- 
ing the ambient photon field. 

In a fluid with r± < 1, bulk Compton drag is most effective at 
the outer scale Lo of the turbulent spectrum. The drag timescale 
is 

Here SV e /c <~ SB/B is the velocity with which the electrons are 
advected past the ambient photon gas by the MHD wave mo- 
tions; and U 1 is the energy density in the photons. Normalizing 
eq. (1 10) to the flow time in the bulk frame, and re-expressing 
the magnetization parameter of the pairs in terms of the bulk 
frame compactness l B and scattering depth r±, 



ct± = — ~3-10, 



one has 



rag 



l 



(111) 



(112) 



r/T c c 

We see that fd rag is comparable to the flow time if t± <~ 1 and 
most of the energy in the magnetic field has already been trans- 
ferred to the photons. 

The timescale for 3-mode wave interactions can, however, be 
significantly shorter than the flow time. In the remainder of this 
section, we examine in more detail how energy is transferred to 
the electrons and positrons through the formation of a spectrum 
of high-frequency waves, and how the heated particles cool. 

6.1. Electrostatic Acceleration vs. Resonant Cyclotron 
Heating of Particles 

A turbulent spectrum of torsional MHD waves is created 
through non-linear couplings between the waves on the outer 
scale Lq. Both compressive (fast) and torsional (Alfvenic) 
waves are created in the fluid, but the torsional waves have 
stronger couplings with each other than do fast modes of a sim- 
ilar frequency (Troischt & Thompson 2004). Alfven turbulence 
has a strong tendency to anisotropy, the wavepackets becoming 
elongated along the magnetic field at high frequencies (Higdon 
1984; Goldreich & Sridhar 1995). This property is preserved in 
the relativistic regime (Thompson & Blaes 1998; Cho 1995). 

The energy that cascades from the outer scale Lo = eo(r/T c ) 
(eq. [109]) can be written as 

B 2 

Ecus — ^caSo (113) 

per unit volume. The coefficient e cas depends on the strength of 
the turbulence. It takes the maximum value e cas <~ £q' (SBq /B) 2 
when the coupling between colliding wavepackets is strong on 
the timescale ~ Lq/c. The three-mode coupling is strong when 

^~^~1 ^ (114) 
SB k\\ B 

8 Measured in the bulk frame of the outflow. 



(Goldreich & Sridhar 1995). Here fcjj 1 is the the length of a 

wavepacket parallel to the background magnetic field, and k~^ 
is the transverse size. We can then set a lower bound on this 
coupling parameter from the requirement that a cascade have 
time to develop on the (bulk frame) timescale r/T c . When the 
parameter is small, the effects of subsequent collisions will ac- 
cumulate as a random walk. A wavepacket of length <~ Lo accu- 
mulates a distortion [A(SB)/5B] 2 ~ SQ l (k±5B /k\\B)l in a time 
r/T c c. This distortion is of the order of unity if (k±5B/k\\B)o > 

e l J 2 . Taking (k±/k^) ~ (V rec /c)~ l , and e from eq. (109), one 
sees that a significant fraction of the total magnetic field energy 
will transfer to high-frequency torsional waves if 

' SB A\ „ . . ,~-3 r c2 {Vrjc\ ( A-r B /c ' 
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This inequality is satisfied if 5B /B > 0.2. 

If coupling parameter (1 14) is smaller than unity at the forc- 
ing scale, it grows with wavenumber and the cascade quickly 
enters the 'critically-balanced' regime where k±SB/knB <~ 1 
(Goldreich & Sridhar 1995; Cho 2005). One then is led to a 
simple result that has considerable importance for the radiative 
emission from the turbulent medium (Thompson & Blaes 1998; 
Lyutikov & Thompson 2005): the turbulent spectrum must be 
cut off at a frequency that lies well below the gyrofrequency of 
the electrons and positrons. 

A torsional wave transports charge along the background 
magnetic field at a rate 5 J = ck±5B/4ir. Balancing this against 
the maximum conduction current that can be supplied by the 
ambient electrons and positrons, one has 

k ± SB hB 

(U6) 

At higher current densities, a displacement current develops 
through the Maxwell equation 
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Only a small imbalance on the right hand side is sufficient to 
generate an accelerating electric field that deposits energy in 
the particles. 

The charge-starvation scale (1 16) is therefore 
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Referenced to the outer scale Lo, this is 

h *- U ~aMB QED y 

Here t± = o- T n±Lo is the scattering depth of the pairs across 
a distance Lo, ctem = 1/137 is the fine structure constant, and 
Bqed = m e c 3 /eh = 4.41 x 10 13 G is the QED magnetic field. 
The field in the dissipation zone 8 is weaker than Bqed by more 
than 10 orders of magnitude (eq. [74]). One therefore has 
^|| starve ~ 10 12 t± Lq 1 . One sees that torsional wave turbulence 
can be supported in a GRB outflow over a very wide range of 
frequencies. 

In general, we will be considering the case where the rest en- 
ergy density in electrons and positrons smaller than the energy 
density of the background magnetic field, 
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In this regime, the velocity of a torsional wave (the Alfven' 
velocity Va) is generally close to the speed of light. However, 
Va is significantly smaller than c when the wave is so strongly 
sheared that k± > uop e /c. The cold plasma dispersion relation 
reads 

uj~V A k h (121) 



at low frequencies uo <C up e 
speed is 



(47rn±e 2 /m c ) 1 / 2 . The Alfven 
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including the effects of particle inertia and a strong gradient 
perpendicular to the background magnetic field. Here cr -1 = 
(j ~± + cr" 1 is the total magnetization parameter, where 
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refers to the ions. (The derivation of equation [122] is similar to 
that given by Arons & Barnard 1986 for sheared Alfven waves 
in the infinite-B limit.) 

In this regime, the energy in the wave is divided nearly 
equally between the transverse fluctuation in the electromag- 
netic field, and a longitudinal oscillation of the electrons and 
positrons. Namely, 
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At a perpendicular wavenumber k± > V e /u>p e , the phase speed 
of the mode along the magnetic field drops below the speed of 
light, and the mode is Landau damped on the longitudinal mo- 
tion of the electrons. 

When the ions are warm, the excited modes are electron- 
supported at perpendicular wavenumbers larger than the inverse 
ion gyroscale, k± ^> eB/m p cV p . It should be emphasized that 
the dispersion relation of these electron-supported modes is not 
the usual whistler dispersion relation if the kinetic energy den- 
sity of the ions is smaller than the magnetic energy density, 
jnijniV 2 <C B 2 /Stt. The dispersion relation is instead Alfven- 
like, and differs from (121), (122) only in the absence of the 
factor (1 + <T~ 1 )~ 1//2 . Another difference with non-relativistic flu- 
ids is that modes of both (circular) polarizations are supported, 
and so we conjecture that the mode-mode interactions are sim- 
ilar to those of lower-frequency magnetohydrodynamic waves. 
(More generally, the electron-supported modes will have an 
Alfven-like dispersion relation over some range of wavenum- 
bers k± <ujp e /c as long as n e m e c 2 -C B 2 /ir. Our considerations 
are therefore relevant to a non-relativistic medium in which pc 2 
exceeds B 2 /8ir, but not by a factor exceeding <~ m p /m e .) 

We will generally be interested in the case where there is a 
constant energy flux through a critically-balanced spectrum of 
torsional waves, 
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If the coupling parameter (1 14) is unity at the outer scale k^ <~ 
Lq\ one has k ± L Q = ((5B /5)" 1 (/tjjL ) 3/2 (Goldreich & Sridhar 



1995). The scale k± ~ uip e /c is therefore reached at a parallel 
wavelength 
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Note that the dependence of Va on the total magnetization pa- 
rameter a cancels here. Equation (127) can be re-written as 



k\\ shear U) ' 



Uo) 



2/3 



B 



QED 



ir) ,128) 



using eq. (119). This transition occurs at a wavenumber 

k\\ shear ~ 10 fc|| starve (129) 

when t± <~ 1, er± <~ 10, SB Q /B <~ 1, and the magnetic field is 
given by eq. (74). 

A torsional wave of frequency uj will resonate with the gyro- 
motion of an electron (Lorentz factor j e and speed /3 e || parallel 
to the background magnetic field) only if 

eB 

Je(u-f3 e \\ck\\ stmYf ,) = ± . (130) 

m e c 

We have found that V A <C c at such high frequencies, and so the 
resonant energy is 

le\Pe\\ I ~2a±. (131) 

If all the charges obtained this Lorentz factor and did not lose 
energy to radiation, then the kinetic energy density in particles 
would be comparable to the energy in the background magnetic 
field, 

{ ( 7 2 -lf 2 )m e c 2 n ± ~^. (132) 



This shows that resonant heating of the particles will be pre- 
vented if one of two conditions are satisfied: 

1 . The particles cool significantly in the time required for 
the MHD wave energy to cascade to small scales; 

2. The net wave energy deposited in particles is small com- 
pared with the background field energy, (SB) 2 <C B 2 . 

6.2. Heating of Ions by Alfven Waves 

It is possible in principle for resonant heating of the ions 
to limit the range of wavenumbers over which left-handed tor- 
sional waves can be supported. This effect can be neglected if 
the wave frequency at the scale (118) is smaller than the ion 
gyrofrequency, 
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We take the ions to be protons for simplicity. In fact some sig- 
nificant admixture of heavier ions with a charge/mass ratio ~ j 
is expected in most circumstances. This additional component 
is mainly helium if the baryons are derived from neutronized 
material in the central engine (e.g. Beloborodov 2003); but 
could be carbon and oxygen if there is mixing between a rel- 
ativistic jet and the outer layers of a Wolf-Rayet star. 

To evaluate this expression, it suffices to note that Va oc k~^ 
when k± > up e /c (eq. [122]). We assume that the cascade 
remains critically balanced (i.e., set the coupling parameter 
k±5B/k\\B <~ 1), and that the Kolmogorov energy flux (126) 
is conserved. Then we find SB ~ constant and k^ oc k± when 
k± > uope/c. The left-hand side of eq. (133) is therefore equal 
to its value at the transition scale kT, l h . Defining 
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we have 
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This expression simplifies in a non-relativistic medium where 
pairs are absent. After expressing the free-electron density 
in terms of the magnetization parameter a e = B 2 /Snn e m e c 2 = 
(m p /2m e )(VA/c) 2 , one has 
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We reach the interesting conclusion that the electrons and 
positrons will experience non-resonant, electrostatic heating 
even in a fluid where pc 2 > B 2 /%tt (but not by a factor exceed- 
ing ~ m p /m e ). Non-resonant heating occurs generically if the 
magnetic field carries even a modest fraction of the outflow lu- 
minosity. 

6.2.1. Implications for Black Hole Coronae 

We now consider the implications of these results for the 
damping of magnetohydrodynamic turbulence in black hole 
coronae. An additional damping mechanism, Compton drag 
by the ambient radiation field, is available to MHD waves in 
a black hole accretion disk. We argue that, in this case, the 
dominant damping mechanism remains the turbulent cascade. 
For the turbulent layer at the base of the black hole corona, 
we take V A /c - 0.1 and B~7x 10 6 (M BH /10 Mg)" 1 / 2 G at 
r <~ 10GMbh/c 2 (the magnetic field which supplies the accre- 
tion torque: eq. 5.9.10, Novikov & Thorne 1973). This gives 
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Thus, electrostatic heating becomes effective before resonant 
absorption on the ions (see also Thompson & Blaes 1998). (A 
relatively small proportion of the cascade energy can be trans- 
ferred to the ions by Landau damping on the gyrational motion 
of the ions 9 at an intermediate wavenumber k± <~ eB/m p cV p : 
Quataert & Gruzinov 1999). 

Now let us examine the Compton drag of the X-ray (and 
UV) radiation field acting on the bulk magnetohydrodynamic 
motions. The drag timescale is fd rag ~ r im p cj\uiU 1 in a non- 
relativistic medium whose inertia is dominated by the ions. 
Comparing fdrag with the non-linear damping timescale f NL = 
l/k\\V A gives 

fdrag _ 3(c/V A ) ( B 2 



(138) 
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where we have defined t t (A||) = n e a-Yk~^ . It has been sug- 
gested that the X-ray photons emerging from an optically 
thick, radiation-dominated disk could be upscattered into a non- 
thermal tail by these bulk Alfven motions even close to the disk 
mid-plane, where the scattering depth tt ^> 1 (Socrates, Davis, 
& Blaes 2004). Setting aside the question of whether the radia- 
tion density will be as high as the thermally unstable Shakura- 
Sunyaev disk solution would imply, we note that the formation 
of a non-thermal X-ray continuum requires a significant frac- 
tion of the energy flux to be carried by the magnetic field. That 
is, one requires 

B 2 c 

ttV a > —U 1 (139) 
87T 3t t 

9 I thank Yoram Lithwick for a discussion of this point. 

10 Of photons with rest frame energies £ 7 < m e c 2 /-y e . 



at a vertical Thomson depth r T > 1 . Substituting this expression 
into eq. (138) gives 
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To evaluate this expression, note that 
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where h is the vertical scale height of the dissipating layer. This 
shows that bulk Compton drag is likely to be less effective than 
3-wave couplings at damping MHD turbulence at the base of 
a black hole corona. We emphasize that the net effect is still 
to deposit the wave energy directly in the electrons, but by the 
mechanism of electrostatic heating. 

6.3. Compton Cooling 

We showed in the preceding section that electrostatic heat- 
ing of the light charges in a turbulent magnetoplasma can occur 
gradually, through many small impulses. Heating can, there- 
fore, be compensated by an incoherent cooling process. In 
this section, we compare the relative importance of inverse- 
Compton and synchrotron cooling. 

Compton drag will damp the motion of the charges on a 
timescale 

3m e c ( 3 \ r 
hc = - —= )w- (142) 
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where t/ 7 is the ambient radiation energy density and 
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is the corresponding compactness in the bulk frame. 10 The en- 
ergy equation 
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has the equilibrium solution 
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The corresponding Compton parameter can be expressed in 
terms of the Thomson depth t± = n±a-\r/Y c as 



yc = Mft)T±=e a 



8tt?7^ 
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This parameter can be large early on in the process of turbulent 
damping, but must quickly saturate at a value yc < 1, which a 
significant fraction of the available magnetic energy density has 
been transferred to the radiation field (Thompson 1994). The 
heated particles can sustain relativistic energies only if r± < 1 . 

Modulations in the volumetric heating rate will occur over 
a range of lengthscales and timescales in the fluid. The above 
estimate of the equilibrium Lorentz factor does not apply when 
the cascade timescale becomes shorter than the inverse Comp- 
ton time-scale. We examine this distinct regime of 'flash' heat- 
ing in §6.4. 
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Let us now estimate the relative importance of synchrotron 
cooling. When t± > 1, the electrons and positrons are 
maintained at a sub-relativistic temperature. Direct cyclo- 
synchrotron emission is suppressed by the large optical depth 
at the cyclotron resonance, 

7T 2 en e (r/T c ) t± 

Tcyc ^ , r. (147) 
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Even when the torsional wave spectrum is cut off at a fre- 
quency much lower than the electron gyro frequency eB/m e cj e , 
the light charges will still acquire momentum perpendicular to 
the magnetic field by upscattering ambient photons. An elec- 
tron moving relativistically along the magnetic field receives a 
transverse momentum 

E' 

p_L=sin6>'-^ (148) 

c 

when it scatters a photon into a direction 9' measured with re- 
spect to the field. Here ~ 7 e £' 7 ( 1 - cos 9) < m e c 2 is the pho- 
ton energy in the electron rest frame, and 9 its angle of propa- 
gation in the fluid frame. The perpendicular momentum gained 
per scattering is obtained by averaging over 9 and 9' , using the 
differential cross section daj/dfl' = (3ot/167t)(1 + cos 2 6*') 
and assuming an isotropic photon distribution in the fluid frame. 
One finds 
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at fixed energy E 1 . The effects of successive scatterings accu- 
mulate as a random walk, giving 
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over a timescale r/T c c. 

The synchrotron cooling rate can now be calculated as 
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Comparing with the inverse-Compton power gives 
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In this expression the radiative compactness t~ (eq. 143) 
has been rescaled to the total magnetic energy density, Ib = 
£-y(B 2 /8TrU-y). The inverse-Compton emission from an ensem- 
ble of charges of energy j e is beamed into an angle 

(Pl) l/2 
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The anisotropy of the electron distribution, and the relative 
output in synchrotron and inverse-Compton radiation, depends 
on the spectral distribution of the ambient radiation field. When 
the radiation field has a Wien distribution with a bulk-frame 
temperature T r , then 
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If the observed temperature T r0 bs is typical of the spectral peak 
of gamma-ray bursts, ^B^obs = \T c k^T r < m e c 2 , then the bulk- 
frame temperature will be far below m e c 2 /k B . In this case, the 
equilibrium pitch angle of the radiating electrons is small, 
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The synchrotron power is also small relative to the inverse- 
Compton power as long as l B < T c (which is indeed the case 
when the deceleration is tied to the end of pre-acceleration; eq. 
[46]). 

Similar conclusions are reached when the target radiation 
field has a high-energy tail. The scattering is in the Thomson 
regime as long as £ 7 < m e c 2 /-f e . At higher energies the cross 
section is suppressed by a factor ~ ^(m e c 2 /E' )ln(E' /m e c 2 ). 
When the spectrum above the observed peak energy E p&l ± is a 
power law with energy index (3 = -1 up to a maximum energy 
£max > r c (m e c 2 /j e ), the average in eq. (150) becomes 

(E 2 ) m e c 2 l + (3/8)ln 2 ( 7c £ max /r c m e c 2 ) m e c 2 



where the last equality holds for 7, 
electrons saturates at 
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1 . The pitch angle of the 
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and the ratio of synchrotron power to inverse-Compton power 
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The above expressions for the r.m.s. pitch angle apply to 
a charge that is heated continously over the expansion time 
~ r/T c c. The pitch angle is smaller if a charge is accelerated 
suddenly along the magnetic field, and then cools passively on 
a timescale ~ (^ 7 7 e ) _1 r/T c c. Since {p\} grows linearly with 
time, one has 



^synch (p±) 



Ejc " £ 7 hem e c) 2 ' £^ 2 lnT c (/? 1} ' (159) 

The r.m.s. pitch angle of the cooling particle never grows larger 
than 1 / 7e, and the emission cone of an ensemble of particles is 
comparable in width to that of a single particle. 

These expressions are easily generalized to softer gamma-ray 
spectra, with an energy index -2 < (3 < -1. One finds 
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in the case of continuous heating; and 
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in the case of sudden heating followed by passive Compton 
cooling. 

6.4. Flash Electron/Positron Heating 

We have argued that the relativistic outflow is likely to con- 
tain a variable magnetic field, which switches sign on a length- 
scale that is small compared with cAf (but still much larger 
than the size of the neutronized torus that feeds the central black 
hole). As a result, it is plausible that stochastic bursts of dissi- 
pation (driven by magnetic reconnection) occur on timescales 
that are small compared with At (that is, on lengthscales that 
are small compared with r/T c in the frame of the contact). This 
means that a fraction (SBq/B) 2 of the energy of the shock rel- 
ativistic outflow can be transferred locally to the electrons and 
positrons quickly enough that they reach the equilibrium energy 

' SBci n 2 
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before they are able to cool. This sets an upper bound on 
the size Lq of the heated region. One requires that Lq/c < 
*cooi((7e)) = 3m e c/4(j e }ajU-y, or equivalently that 

Lq ^ 1 t± /SBqY 2 

< — — = — I — - I . (163) 
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Here we have related the magnetization parameter of the pairs 
to the scattering optical depth via <t±= I b /t±. 

7. BEAMING OF THE INVERSE-COMPTON RADIATION IN A 
TURBULENT MAGNETOFLUID 

We have found that a relativistic shell of ejecta begins to dis- 
sipate strongly when the compactness of the radiation stream- 
ing across the forward shock has dropped below a characteristic 
value (§3.1). This results in a characteristic radiative compact- 
ness I' within the moving shell (eq. [46]), and a characteris- 
tic optical depth to scattering between the reverse and forward 
shocks (the sum of eqs. [54] and [177]). Because the dissi- 
pation is concentrated in a narrow range of radius, one obtains 
an explanation for the weak evolution of the temporal power 
spectrum that is typically seen within a gamma-ray burst. 

If the gamma-ray emission is triggered by the interaction 
with the external medium, then the observation of many well- 
separated subpulses in some gamma-ray bursts forces a signif- 
icant constraint on the the emission mechanism: the gamma- 
rays must be beamed in the bulk frame (in the frame of the 
contact discontinuity). By contrast, high-frequency variability 
from internal shocks is smoothed out due to the curvature of the 
ejecta shell, unless the internal shocks are occurring well inside 
the radius where the reverse shock passes through the ejecta 
shell (Sari & Piran 1997). Emission on a timescale St <C r/2cT 2 c 
must be localized to a small fraction ~ (2Y 2 c5t / r) 2 of the sur- 
face area of the shell that is visible to the observer. The emis- 
sion occurs at a radius r <~ 2T 2 cAt when it is triggered by the in- 
teraction with the external medium. This means that ~ (At /St) 2 
independent regions of the shell contribute to each time interval 
St, and only a fraction of them must be visible to the observer 
in a highly variable burst. 

What is the cause of this beaming? Reconnection in a high-cr 
magnetofluid can, in principle, create bulk motions as fast as 
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where V 2 = (B'f/Airp 1 (Blackman & Field 1994) and cr" 1 = 
<j~ l + a± . So one could ascribe the gamma-ray variability to 
stochastic variations in the direction of bulk motion (Lyutikov 
& Blandford 2003). It is not clear, however, that significant 
elements of the fluid will be able to accelerate to relativistic 
speeds, as they evidently do in the outer magnetosphere of a 
Soft Gamma Repeater during a giant flare (e.g. Thompson & 
Duncan 2001). In contrast with an isolated neutron star, whose 
magnetic dipole field pressure drops off rapidly with radius (as 
<~ r~ 6 ), one expects the pressure to equilibrate rapidly after the 
relativistic outflow has passed through the reverse shock. An 
element of magnetofluid may therefore feel a strong drag force 
off the ambient fluid as it attempts to move relativistically. 

A quite different beaming mechanism is afforded by our 
model of decaying Alfvenic turbulence (§6). The energy of 
the waves is extracted by electrostatically accelerating electrons 
and positrons along the magnetic field. As a result, the radia- 
tion of these charges is collimated along the local direction of 
the field lines (Fig. 4). 



The beaming pattern depends in a subtle way on how the 
charges are heated. The charges start out with small per- 
pendicular motion in the case of flash heating (where energy 
is transferred from the torsional waves on a timescale much 
shorter than the cooling time; §6.4). Their radiation is there- 
fore beamed into a solid angle <~ tt/j 2 . Repeated inverse- 
Compton scatterings impart some gyrational motion. By the 
time a charge has lost half of its initial kinetic energy, the pitch 
angle has increased to (p 2 L ) 1 ^ 2 /j e m e c ~ l/^ e (lnT c )^ 2 (in the 
case of hard high-energy photon spectrum; eq. [157]). Here 
7 e ~ <j±(5B[ ) /B') 2 depends on the fraction of the magnetic en- 
ergy that is transferred quickly to the charges (eq. [162]). The 
beaming pattern of a large collection of charges does not change 
significantly from that of a single charge. 




FIG . 4 . — Torsional MHD waves in a magnetofluid transfer energy to the light 
charges by forming a high-frequency spectrum of current fluctuations, and an 
accelerating electric field. The charges cool primarily by inverse-Compton 
scattering ambient photons along the local direction of the magnetic field. The 
observed duration of a gamma-ray spike decreases with energy (eqs. [170], 
[169]). The scaling with energy depends on whether cooling is faster than the 
shift in beaming direction that is caused by the long-wavelength undulations 
of the magnetic field. 

The heating may, alternatively, be slow enough to be bal- 
anced by inverse-Compton drag. The electrons and positrons 
will then maintain a pitch angle (157) - but at a somewhat 
smaller Lorentz factor, (7 2 ) 1 / 2 ~ (3o-±/4£' 7 )V 2 (eq. [145]). 
The radiation of many charges is now beamed into a larger 
solid angle tt9' 2 ~ ir (p 2 ± ) / \^ e m e c) 2 <~ 7r^' /7 e lnr c . One sees 
that beamed radiation is a direct consequence of the damping 
of Alfvenic turbulence. 

The beaming pattern of the gamma-rays, as seen by the ob- 
server, is sensitive to the presence of slow undulations of the 
contact, through an angle larger than the beaming angle 59' . 
The inverse-Compton radiation is beamed in the directions par- 
allel and anti-parallel to the magnetic field. An observer sit- 
uated in the rest frame of the central engine will see the seed 
thermal photons emitted within an angle ~ 1 /T c of the radial 
direction. The inverse-Compton photons will, by contrast, be 
visible only in two spots of angle diameter <~ 69' /T c that are 
displaced by an angle ~ 1/2T C from the center of the visible 
patch of the ejecta shell. 

7.1. Correlation between Variability Timescale and 
Photon Energy 

We now consider the temporal pattern of the beamed radi- 
ation that is emitted by the heated pairs. A torsional wave 
wave of an amplitude £ and a period P' = 2n/(kuc) imparts a tilt 
SB' /B' = k\\ £ to the magnetic field. This tilt varies with time, 
and so the wave motion forces the radiation beam to to sweep 
past an observer on the timescale 

59' 59' P' 1/2 

StB ~m^JdF = 2«(6B>/B>) Cck \\ (165) 
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(in the bulk frame). The fluctuation timescale is shorter in the 
observer's frame, 5t B = 5t' B /T c . 

This expression for 5t' B holds when SB' /B' > 59', which is 
possible only over a certain range of wavenumbers in a criti- 
cally balanced Alfvenic cascade. The conservation of energy 

—1/2 

flux (eq. [126]) implies that SB' /B' oc ' . The minimum 
fluctuation timescale is 

aT^U^J' (166) 

given a wave amplitude SB' at an outer scale Lq = so (T e cAf ) in 
the bulk frame. 

Two other timescales for variability are relevant here when 
the particles undergo flash heating. The motion of the charges 
along the magnetic field is damped on the timescale 

(as seen by an observer who is at rest with respect to the central 
engine). Finally, the beaming angle 89' broadens in response 
to repeated Compton scatterings. In particular, the r.m.s. pitch 
angle becomes comparable to ~ 1 on a timescale 

aT^- (168) 

(This expression assumes a hard photon spectrum, (3 = -1, and 
is easily generalized to softer spectra using eq. [160].) 

All of these timescales are potentially resolvable in a gamma- 
ray burst light curve. By contrast, if the emission mechanism 
were synchrotron radiation, then the cooling time would be far 
too short to agree with the observed width of the gamma-ray 
sub-pulses. One requires a Lorentz factor j e > 10 4 to gen- 
erate photons with a observed energy of 100 keV or greater 
(eq. [74]). The corresponding synchrotron cooling time is 
At/£'^ e <lO- 4 At. 

Each of these variability timescales become shorter as the en- 
ergy of the inverse-Compton photon increases. When photons 
from the advected thermal peak are the principal seeds, one has 
£, 7^7 c 2£ 'peakand 



8t coo i, Stg oc E ' 1/2 ; 



(169) 



St B oc 59' oc E~ X I A 



(continuous heating); 
St B oc 59' oc £" 1/2 (flash heating). (170) 

Which of these scalings is most applicable depends on the speed 
of Compton cooling relative to the change in beaming direction 
that is driven by long-wavelength undulations of the magnetic 
field. It is interesting to recall, in this regard, that some gamma- 
ray bursts are composed of spikes which show noticeable asym- 
metry between rise and decay - the so-called FRED behavior 
(e.g. Fenimore et al. 1996). Other bursts are highly variable 
but composed of spikes with no clear asymmetry between rise 
and decay. We associate this second type of variability with 
fluctuations in beaming direction on a timescale that is shorter 
than the cooling time. FRED pulses demonstrate the scaling 

between width and frequency that is characteristic of cooling, 

—l li 

5t oc E 7 (eq. [170]; e.g. Fenimore et al. 1995; Norris et 
al. 1996). An analysis which separates out the two types of 
variability has not yet been done, but would be illuminating. 

' 1 We ignore spectral evolution for the purpose of this comparison. 



7.2. Pair Creation: Effects of Beamed Emission and 
Inertia of the Breakout Shell 

The ejecta are photon rich during the deceleration process. 
The number of photons emitted per unit area inside radius r is 

This can be normalized to the electron column by 1 1 
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Substituting eqs. (30), (34), and (36) gives 
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The shocked outflow will, itself, have a modest optical depth 
to scattering (eq. [54]). The larger contribution to the optical 
depth t± comes from gamma-ray collisions within the shell. An 
estimate can be obtained directly from the observed gamma-ray 
spectrum, under the assumption that the high-energy power-law 
extends up to the pair-creation threshold <~ m e c 2 in the emitting 
frame (e.g. Baring & Harding 1997). The standard estimate as- 
sumes isotropic emission of the gamma-rays in the bulk frame. 

Given an observed peak energy £p ea k> the density of photons 
of energy m e c 2 in the bulk frame is 

'Y c m e c 2 ^ 



n'Jm e c 2 ) < 



47rr 2 r e £ peak cln(OT e c 2 /£p eak ) \ E ] 



^peak 



(174) 

We specialize here to a hard power-law with a photon energy in- 
dex -1 extending from energy £ pea k up to a bulk frame energy 
~ m e c 2 . The density (174) can be re-expressed as 



a T n'(m e c 2 ) - 



(175) 



ln(m e c 2 /£ P eak) ' 

where £' y is the bulk frame compactness of the photons, defined 
analogously to eq.(43), 



l' = 



(176) 



^m e Cr 

The density of created pairs is n'± ~ 2e 77 ctt [« 7 (m e c 2 )] 2 (r/T c ), 
where e 77 normalizes the frequency-averaged photon collision 
cross section to err, and is e 77 ~ 0.1 for a (3 = -1 photon energy 
index (Svensson 1987). The scattering depth in the region of 
shocked relativistic material behind the contact is then 



T±(r<R c ) = e~, 



2(0 2 



In (r c m £ ,c 2 /£ P eak) 
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Here we have substituted eq. (45) using e+ = 0.25 at the radius 
where half the ejecta are shocked. 

There is, however, considerable reason to doubt this assump- 
tion of isotropic emission. Some bursts show strong variabil- 
ity, which in the framework advanced in this paper demands 
that this emission be beamed in the bulk frame (§7). In such a 
situation, the production rate of pairs rises considerably. Sup- 
pose, for example, that the emission is collimated within cones 
of angular width 259' ~2/j e that run anti-parallel (in the bulk 
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frame). We fix the net emissivity per unit volume. The inten- 
sity of gamma-ray photons in the oncoming beam increases by 
a factor 2-k/-k{89') 2 ~ 2^. In addition, the threshold condition 
for pair creation by the two colliding photons, 

£ 7l £ 72 (l-cos6> 12 )>2(m c c 2 ) 2 , (178) 

is more easily satisfied. If the photons are moving anti-parallel, 
#12 ~ 7r, then the threshold energy for pair creation is reduced 
by a factor ~ 1/V2 compared with the r.m.s. value in an 
isotropic photon gas. Assuming the same spectrum as in eq. 
(174), the pair creation rate therefore increases by a factor ~ 2. 

The breakout shell also has a strong effect on the compact- 
ness of the shocked fluid. Inward from the point where T c 
reaches its peak value, the compactness rises rapidly toward 
small r, £'(r) oc r~ 2 (eq. [72]). The proportion of the relativistic 
ejecta which has been shocked decreases only slowly inward, 
as Ar oc r 1 / 3 . The optical depth in pairs is proportional to (£') 2 . 
Combining these scalings, one sees that t± is an extremely 
strong function of the thickness of the shocked relativistic shell, 

t± oc (Ar) -12 . (179) 

In addition, the relativistic jet material can develop a mildly rel- 
ativistic motion with respect to the fragmenting breakout shell. 
This reduces the energy threshold for pair creation by allowing 
some of the gamma-ray photons to side-scatter off the denser 
shell material. 

We conclude that a substantial fraction of the relativistic out- 
flow will dissipate at optical depths r± > 1, when the effects of 
beaming and the breakout shell are taken into account. 

8. IMPLICATIONS FOR THE SPECTRA OF 
GAMMA-RAY BURSTS 

We have explored in some detail the physics of a relativis- 
tic outflow, composed of a gas of seed thermal photons and a 
stochastically reversing magnetic field, which is launched into 
the dense wind of a Wolf-Rayet star. In this section, we out- 
line the implications of our results for the non-thermal spec- 
tra of gamma-ray bursts. We first address the origin of the 
high-energy non-thermal continuum that is a defining feature of 
GRBs. Second, we show how the observations of gamma-rays 
at energies > 1 2 times the peak energy provide a valuable diag- 
nostic of the physics of pair-loading, the density of the ambient 
medium, and the compactness of the zone in which the ejecta 
decelerate. Third, we address the spectral transition that occurs 
between the end of the prompt GRB phase and the beginning 
of the afterglow emission (a good example being GRB 980923: 
Giblin et al. 1999). We then explain how prompt optical syn- 
chrotron emission may be suppressed (without self-absorption) 
in a turbulent magnetoplasma. Finally, we connect our model 
with a previous suggestion (Thompson 1994) that thermal ra- 
diation generated close to the central engine is the dominant 
source of inverse-Compton seeds for the non-thermal radiation 
of GRBs. The physical mechanisms that we have examined in 
this paper can also be applied to spectrally harder classes of 
bursts - the short GRBs (Kouveliotou et al. 1993) and the giant 
flares of the SGRs (e.g. Hurley et al. 2005). In these types of 
bursts, one requires that the outflow go into (nearly) free expan- 
sion much closer to the engine than it does in the Wolf-Rayet/jet 
model. 

Two general lessons emerge here. First, a close correspon- 
dence between the peak of the seed thermal spectrum and the 
peak of the final GRB spectrum arises much more easily if the 



heating of the pairs is spatially and temporally distributed, in- 
stead of being localized at a strong shock wave (or a few such 
shock waves within each causal patch of the outflow). Even in 
a situation where the magnetic energy density greatly exceeds 
the rest energy density of the radiating particles, rapid Compton 
cooling allows the mean particle Lorentz factor to remain close 
to m e c 2 in the fluid rest frame. 

Second, single-box models of the radiation spectrum can pro- 
vide some subtle insights into the interplay between different 
radiative processes (e.g. Pe'er & Waxman 2004b), but ulti- 
mately the GRB emission problem is one of radiative transport. 
The prompt GRB emission continues only as long as the ther- 
mal seed photons continue to overlap with the dissipating mag- 
netized shell. The prompt 100 MeV component of the GRB 
spectrum is a byproduct of the dissipative processes occurring 
two distinct regions: the inverse-Compton emission behind the 
contact discontinuity; and the side-scattering of these gamma- 
rays in the Wolf-Rayet wind. The fluid that passes through the 
reverse shock builds up pairs by collisions between gamma- 
rays, but this pair creation is extended in time. As a result, 
the mean energy per particle has a strong negative gradient 
away from the reverse shock. A high-energy spectral tail is the 
byproduct of this inhomogeneous distribution of particle ener- 
gies. 

8.1. Power-Law Spectral Tails above the Peak Frequency 

The spectra of gamma-ray bursts are characteristically non- 
thermal, and in many cases appear to have a power-law shape 
above the peak energy £p e ak- BATSE burst data do not charac- 
teristically cover a wide energy range above the peak energy; 
but some Ginga bursts with lower peak energies have high- 
energy power-law tails that cover more than two orders of mag- 
nitude in frequency (Strohmayer et al. 1998). 

Compton upscattering of the seed thermal photons will pre- 
serve the relation between peak energy and isotropic luminosity 
if i) the optical depth through the dissipating shell is r± > 1 ; and 
ii) the particle heating is distributed broadly through the fluid. 
We have seen that the thermal photons follow the Amati et al. 
(2002) relation - both in slope and normalization - if they are 
generated in the turbulent core of the jet near the photosphere 
of the Wolf-Rayet star (eq. [8]). And the scattering depth can 
attain values near unity in an outflow of isotropic luminosity 
10 51 L re i5i ergs s" 1 and duration lOAfi s, moving in a Wolf- 
Rayet wind with mass-loss rate 10~ 5 M w -^Mq yr -1 (eq. [177]). 
This requires a hard photon spectrum (energy index (3 = -1), 
which is consistent with a significant majority of BATSE bursts 
(Band et al. 1993) and essentially all of the sample of Ginga 
bursts analyzed by Strohmayer et al. (1998) (excepting those 
with high peak energies in which the high-energy tail was not 
fully sampled). Even when M w < 10~ 5 M Q yr" 1 , the optical 
depth can be strongly enhanced by the residual inertia of the 
breakout shell, and by angular inhomogeneities in the outflow 
that allow side-scattering of photons off slower, denser material 
(§7.2). 

A net optical depth t± > 1 in the dissipation zone is required 
for several reasons. First, the energy radiated by each charge 
is inversely proportional to t±. The energy of the inverse- 
Compton photons scales as t± at small optical depths, and 
so the peak of the inverse-Compton spectrum lies far above 
the peak of the thermal seed photon spectrum. Second, a hard 
inverse-Compton tail dL-JdE^ oc E^ 1 ^ 2 is generated below the 
peak energy if the base of the non-thermal electron spectrum 
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lies at an energy j e ^> 1 . Third, a significant fraction of seed 
blackbody photons must be scattered, so as to avoid a strong 
and localized thermal bump in the transmitted spectrum . 

A lower bound to the luminosity of the seed photons is ob- 
tained by demanding that the first-order inverse-Compton (IC) 
photons carry a significant fraction of the bolometric output of 
the burst. As the heated particles cool off, and the energy den- 
sity of the IC photons may begin to exceed that of the blackbody 
seeds, in which case the first-order IC photons become the dom- 
inant coolant for the heated pairs. Nonetheless, because £' > 1, 
it is possible for the particles to be energized repeatedly - fol- 
lowed by cooling off the ambient radiation field - so that the 
instantaneous energy remains close to j e <~ 1 . Starting with a 
compactness £' bb in seed photons, it is possible to radiate a to- 
tal energy £' ~ (£ bb ) 2 this way. The seed photons must carry a 
minimal fraction of the outflow energy 

l f > (£'r l/2 (180) 

1/2 

to avoid the appearance of a prominent £y tail in the soft 
gamma-ray spectrum. 

The existence of some bursts (or subcomponents of bursts) 
with very soft high-energy spectra (Pendleton et al. 1997) is, 
of course, fully consistent with this model. The thermal seed 
photons will not be significantly reprocessed if the Wolf-Rayet 
wind has a mass-loss rate much less than ~ 10~ 5 M Q yr" 1 . At 
the same time, one can expect that the optical depth will con- 
tinue to drop as the ejecta expand, so that the later parts of a 
burst will have a tendency to be spectrally softer. So there is 
a general expectation that NHE si/fr-components of bursts will 
be concentrated toward the end of the burst. A good example is 
GRB 920622B (Fig. 2a in Pendleton et al. 1997). 

There is, nonetheless, a basic tension in this model between 
maintaining the £p ea k-£iso relation, and allowing rapid variabil- 
ity in a burst light curve. It should be emphasized that beam- 
ing of the inverse-Compton photons will be maintained when 
t±> 1, because the motion of the electrons and positrons is 
collimated along the direction of the magnetic field. Nonethe- 
less, strong beaming cannot be maintained below the ther- 
mal peak energy at t± > 1, because the Compton parameter 
yc = \t±{01) is limited to a value yc < 1. The relative ab- 
sence of bursts with smooth light curves at the highest measured 
fluxes (Norris et al. 2005) may be connected to the modest de- 

-1/2 

crease in t± that is expected at higher luminosities, t± oc L rel . 
The hypothesis that greater variability corresponds to lower op- 
tical depths is testable: more variable bursts should be spec- 
trally harder, on the average, than the Amati et al. (2002) rela- 
tion would imply. 

A hard gamma-ray spectrum depends on the presence of in- 
homogeneities in the outflow. We now consider three possibili- 
ties that are suggested by this model. 

8.1.1. First-order Fermi Acceleration of Pairs at a 
Mildly Relativistic Reverse Shock, followed by 
Inverse-Compton Cooling 

A non-thermal positron tail extending upward from a min- 
imum energy 7 e <~ 1, and having a significant optical depth 
r ± ^ 1) would generate a hard, high-energy tail to the gamma- 
ray spectrum by cooling off the seed thermal photons (Thomp- 
son 1997). The asymptotic Lorentz factor that the jet material 
attains after escaping the Wolf-Rayet star (eq. [12]) is barely a 
factor 2 larger than the equilibrium Lorentz factor of the contact 



during the deceleration phase (eq. [39]). This means that the re- 
verse shock wave is mildly relativistic, and the mean energy of 
the thermal pairs behind the shock is indeed (l-2)m e c 2 . How- 
ever, our deceleration model suggests that the optical depth of 
the pairs crossing the reverse shock, and the energy they carry, 
are too small to produce the desired effect. 

Collisions between gamma-rays that flow upstream across 
the reverse shock will generate a somewhat smaller depth, 
t± = 0.1r±_i, than they do in the downstream fluid . (This 
optical depth is sensitive to the amount of cooling and com- 
pression that occurs downstream of the shock. The bulk motion 
of the fluid and the beaming of the emission downstream of 
the reverse shock both have the effect of suppressing the pair 
creation rate upstream.) This means that the pairs carry only a 
fraction 

e± ~ y ~ 10- 2 t±_! (181) 

of the outflow energy when they hit the reverse shock. 

The proportion of the outflow energy carried by the ions, 
£b = rfccrit/Tfc, can easily be larger than s±, even though the 
number of positrons exceeds the number of ions by a factor <~ 
(m p /m e )e±/eb- Particle-in-cell simulations of pair-dominated 
shocks show that the shocked positrons can acquire a non- 
thermal tail by absorbing the cyclotron waves emitted by the 
ions (Hoshino et al. 1992). However, the ions are only mildly 
relativistic behind the reverse shock, and so the frequency of 
the ion ring is too low to resonate with the gyromotion of 
the positrons. (While Fermi-acceleration of the ions is possi- 
ble, such a non-thermal ion tail would presumably be a much 
weaker emitter of ion cyclotron waves.) 

Two other problems with this mechanism present them- 
selves. First, rapid variability in the gamma-ray emission is not 
possible at a reverse shock that is triggered by the deceleration 
off the ambient medium (Sari & Piran 1997). And, second, 
non-thermal positrons of an energy "f e <~ 10 2 e^ ^ 2 will emit 
synchrotron photons in the optical band (§8.4). The observed 
upper bound of <~ 10" 4 - 10" 3 to the fraction of the bolometric 
output that is carried in the optical-IR band (Akerlof et al. 1999; 
Vestrand et al. 2005) then requires that e' B < 10 -4 upstream of 
the reverse shock. 

For all of these reasons, we turn to other mechanisms for 
energizing pairs behind the reverse shock, that are special to 
strongly magnetized plasmas. 

8.1.2. Gradient in Electron Energy due to 
Gradual Pair Loading 

The scattering depth builds up after the relativistic outflow 
passes through the reverse shock. The minimal optical depth of 
the entrained electrons and ions is less than (m e /m p )£' < 10" 2 ; 
the optical depth in pairs is roughly an order of magnitude 
larger. 

There is, therefore, a distribution of magnetization parame- 
ters a± in the fluid between the reverse shock and the contact. 
As pairs are created by photon collisions, the shocked mag- 
netofluid will also be heated continuously (by reconnection of 
the reversing magnetic field; §6.4). This leads to a steady de- 
crease in the mean energy (eqs. [145], [162]) of the heated 
electrons and positrons (Fig. 5). 

The energy of the inverse-Compton photons decreases in par- 
allel with (162). To illustrate the effect, we take the pair den- 
sity to incease as a power of time, a± oc B 2 /n± oc t~ a , and 
the wave energy density to decrease as SBl oc r 13 . One then 
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has 7 e oc f" (Q+,3) during repeated episodes of flash heating, fol- 
lowed by Compton cooling. The scaling is slightly different, 
7 e oc j"*"" 1 "/ 3 " 1 " 1 )/ 2 , if the heating is slow enough to be balanced 
continuously by Compton cooling. The energy of the Comp- 
tonized thermal photons decreases with distance behind the re- 
verse shock, as E 7 oc t~ 2{a+,3) (flash heating) and £ 7 oc 
(continuous heating). When the pair-creating photons are sup- 
plied by the emission at smaller radius, one has a ~ 1. The 
fluence F\c in high energy photons is proportional to dB^t, and 
so 

dFlC ocr^ocZC (182) 



d\nE. 



7 



7 ' 



where 7 = ( 1 - /3)/2( 1 + 0) (flash heating) and 7 = ( 1 - /3)/(2+ (J) 
(continuous heating). The same spectrum is obtained, indepen- 
dent of the mode of heating, when the tubulent intensity is in- 
dependent of distance behind the reverse shock (/3 = 0). In this 

— 1 11 

case, the integrated spectrum is hard, dF\c /d \nE^ oc E 1 . One 
also recovers the observed scaling (Fenimore et al. 1995) be- 



tween peak duration and frequency, t oc E 1 
the mode of heating. 
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FIG. 5. — The pair density builds up gradually after the magnetofluid passes 
through the reverse shock wave, due to a finite rate of photon collisions. This 
leads to a gradual decrease in the mean energy of the charges that are heated 
electrostatically by decaying, high-frequency MHD turbulence. A high energy 
gamma-ray spectral tail emerges from the Compton cooling of this extended 
particle energy spectrum (eq. [182]). 



An important feature of this mechanism is that the emission 
of the soft photons lags behind that of the hard photons: the 
softer photons are emitted at later times on average as the pair 
density builds up. The presence of such a hard-soft lag is a 
well established feature of a subclass of gamma-ray bursts with 
smooth and easily separable sub-pulses (Norris et al. 2005, and 
references therein). Its existence therefore provides interesting 
circumstantial evidence for this type of feedback between pair 
creation and distributed heating. 

8.1.3. Turbulent Cells Driven by Reconnection 

The turbulence in the outflow is expected to be strongly inter- 
mittent if the magnetic field reverses sign in a stochastic man- 
ner. The size Lq = eo(r/T c ) of the 'cells' can be related to the 
reversal scale of the field through eq. (109). They are typically 
smaller than the causal scale r/T c , so that the expansion of the 
outflow has a small effect on their internal dynamics. 

A high-energy gamma-ray tail can be generated in cells with 
optical depths At± > 0.1, when the photons are multiply scat- 
tered by the bulk turbulent motions (Thompson 1994, 1997). 
While this mechanism could operate in bursts with smooth light 



curves, it is disfavored in highly variable bursts. Strong beam- 
ing of the inverse-Compton photons requires relativistic bulk 
motion in this case, so that the various orders of Compton scat- 
tering should be well separated in frequency. A smooth power- 
law spectrum is not naturally produced unless the dissipation is 
spread over a wider range of radius. 

8.2. Compton Cooling of Pairs behind the Forward Shock 

A pair-rich medium sweeps past the forward shock dur- 
ing the prompt deceleration phase. Although the number of 
positrons in this flow consistently exceeds the number of ions, 
the proportion of the kinetic energy carried by pairs gradu- 
ally declines from a value A4± — 1 <~ 1 at r ^ /?decei- to <~ 
(^decei+/^decei-)~ 1,/2 ~ 0.2 when the reverse shock has finished 
passing through the ejecta shell (eq. [52]). The total energy 
dissipated at radius r scales as <~ (r/Rfe ce \+j i l 1 ; a fraction <~ | 
of this being dissipated behind the forward shock. 

The kinetic energy of the pairs is readily radiated away be- 
hind the forward shock: the cooling time is a small fraction 
~ 1/^T± of the flow time behind the shock. We have argued 
that the magnetic field is likely to carry a smaller fraction (84) 
of the energy behind the forward shock than do the gamma-ray 
photons. The shocked pairs therefore lose energy by inverse- 
Compton scattering the gamma-ray photons (Thompson 1997; 
Beloborodov 2005; Fan, Zhang, & Wei 2005). 

The spectral distribution of these inverse-Compton photons 
provides a direct probe the dynamics of the ejecta shell, and 
of the coupling between the ions and the lighter charges. First 
let us consider the case where the ions decouple entirely. The 
shocked pairs have a mean energy T± = r e /2r am b, and the char- 
acteristic energy of the inverse-Compton photons is 
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At the end of the prompt emission phase (r = /?decei+)> one finds 

£ic(^decel+) 4 / /?decel+ \ ^ 



ipeak 



= lxl0V^5 1 </- / 5 X / /Ari, 
(184) 

which is Ei C /E bTeak ~ 400L re / 151 M M ,_ 5 for e+ = j. This means 
that the inverse-Compton photons have an energy less than 
~ m e c 2 in the frame of the contact only if £ peak < 50 keV. It 
will be recalled that T± oc T~ 3 , and so 

E K oc T^ 6 . (185) 
The energy of the inverse-Compton photons is a strong func- 
tion of the Lorentz factor of the contact, and a modest increase 
or decrease in the external inertia will push E\c to higher or 
lower energies. 

The total power in Comptonized photons at energies E 1 < 
Eic is 
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Unless the input gamma-ray spectrum is soft at high energies, 
this inverse-Compton tail will not be visible, due to the combi- 
nation Of Small factors (A4± - l)(/?decel+/^decel-)~ 1,/2 - 

A rapid coupling between the gyromotion of the protons 
and positrons could cool the protons. The PIC simulations of 
Hoshino et al. (1992) suggest that the positrons will develop 
a hard, non-thermal tail (number index <f ln/V c +/(iln7 > -2) 
when m p n p > m e n e +. However, resonant excitation of the gy- 
romotion of positrons of energy ~ j e + requires that the pro- 
ton ring emit a significant power in cyclo-synchrotron waves 
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smaller radius. In a rapidly rotating source, this thermalization 
radius can be normalized to the light cylinder, R\,b ~ cP/2ir = 
5 x 10 6 (P/msec) cm; whereas in a Soft Gamma Repeater flare 
it can be normalized to the radius of the star. The blackbody 
temperature is then higher than we deduced for a thermaliza- 
tion radius of = 2 x 10 10 cm (eq. [8]); it works out to 
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at a frequency n(eB /T±m p c) = (m p /m e )(eB/"f e +m p c). Thermal 
positrons (-f e + = T±) behind the shock will be heated only if 
0.31^ > (m p /m e )r£, i.e., only if T± > 20 and n > 6 x 10 3 . 
In the case where the proton ring has a small thermal spread, 
which limits the harmonic to a value n < n max <C m p /m e , 
one requires a secondary mechanism for creating non-thermal 
positrons with an energy j e + > n^ ax (m p /m e )T±. Unless the 
non-thermal component of the spectrum is very hard, its en- 
ergy density is a modest fraction of the thermal positron energy 
density behind the shock (which itself is < 10% of the energy 
density of the protons). For this reason, it is plausible that the 
protons will retain a significant fraction of their kinetic energy 
after thermalizing downstream of the shock. Rapid proton cool- 
ing requires a collective process that persists longer than the 
synchro/Compton cooling time of the absorbing particles. 

It is also useful to contrast this calculation with previous 
estimates, which gives some sense of the subtleties involved. 
The cooling of synchrotron-emitting particles is dominated by 
inverse-Compton scattering of the synchrotron photons if the 
energy density of the cooling particles exceeds the energy den- 
sity of the magnetic field and the synchrotron radiation is suffi- 
ciently soft (e.g. Sari & Esin 2001). It has been suggested that 
the prompt gamma-ray photons will be inverse-Compton scat- 
tered at the forward shock, thereby triggering a pair cascade 
(Fan et al. 2005). This calculation however neglected pair cre- 
ation in the external medium and assumed rapid energy equili- 
bration between ions and electrons, thereby yielding (7^) ~ 10 3 
rather than (7,,) <~ 10 behind the forward shock. Pair creation 
in the Wolf-Rayet wind guarantees at least a ~ 10% conver- 
sion efficiency of bulk kinetic energy to inverse-Compton pho- 
tons at the forward shock. A high-energy particle tail is, of 
course, probably present at the forward shock; but the lower- 
energy part of the inverse-Compton spectrum will be dominated 
by cooling particles that start with an energy ~ (7^). 

It has also been noted that electrons emitting optical syn- 
chrotron photons could upscatter the gamma-ray peak photons 
to GeV-TeV energies (Beloborodov 2005). However, Compton 
cooling of these particles cannot, by itself, explain the low op- 
tical efficiency that is observed in GRBs 990123 and 041219a 
(Akerlof et al. 1999; Vestrand et al. 2005) if the magnetic field 
is strong in the gamma-ray emitting region (§8.4). One infers, 
instead, that the particles must be restricted to a more limited 
range of energy and pitch angle than one would expect from 
standard shock-acceleration. 

Compton GRO detected a very bright burst, GRB 941017, 
whose bolometric output was dominated by a high-energy com- 
ponent with an energy index (3 ~ above 1 - 10 MeV (Gonzalez 
et al. 2003). This high-energy emission probably involved the 
inverse-Compton scattering of the prompt burst photons by en- 
ergetic electrons and positrons. This spectrum was harder than 
eq. (186), and is more consistent with continuously heated par- 
ticles. For this reason, and because of the relatively low normal- 
ization of (186), we favor emission from pairs that are heated 
by a reconnecting magnetic field behind the reverse shock. A 
more detailed discussion of this event is given in §8.6. 

8.3. Harder Classes of Gamma-Ray Bursts 

Even when an extended Wolf-Rayet envelope is absent, one 
still obtains a characteristic temperature for thermalization at 
the base of the outflow. The free expansion phase begins at a 

12 We are most interested in the case of anisotropic emission, with a small pitch angle along a magnetic field with a non-radial orientation in the bulk frame. 

The mean frequency of a synchrotron photon therefore transforms as v = Y c v' from the bulk frame to the rest frame of the engine. 
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The normalization here is well above that relation between £p ea k 
and L; so for long bursts with measured redshifts. 

The same conclusion holds for two recent short gamma-ray 
bursts with tentative redshift measurements: GRB 050509b 
(z ~ 0.2; Bloom et al. 2005), and GRB 050709 (z ~ 0.16; 
Price, Roth, & Fox 2005). In both cases, the spectrum is there- 
fore consistent a binary neutron star merger as the source of 
the burst; and, in the case of GRB 050509b, the energetics is 
marginally consistent with an extragalactic SGR flare. (The 
accretion-induced collapse of a white dwarf following a binary 
white dwarf merger could give rise to the birth of a young, ac- 
tive magnetar in a galaxy with a low star formation rate: King, 
Pringle, & Wickramasinghe2001; Thompson & Duncan 1995.) 

8.4. Suppression of Optical-IR Synchrotron Emission 

The presence of the Wolf-Rayet wind forces the relativistic 
ejecta to decelerate at a radius where electrons of an energy 
7e <~ 10 2 will radiate in the optical band behind the reverse 
shock (eq. [75]). Observations of direct optical-IR emission 
from two gamma-ray bursts (990123: Akerlof et al. 1999; and 
041219a: Vestrand et al. 2005; Blake et al. 2005), and upper 
bounds on other bursts (Akerlof et al. 2000) therefore allow 
stringent constraints to be placed on the particle distribution 
in the outflow. The fraction of the bolometric power released 
in the optical band in these two bursts was ~ 10~ 4 -3 x 10~ 3 
in different parts of GRB 990123, and - 1 - 3 x 10" 4 in GRB 
041219a. Electrons/positrons of an energy j e ~ 10 2 must, 
therefore, contribute only modestly to the energy density in the 
emitting region; or have small pitch angles. 

Taking the simplest case of mono-energetic electrons with 
a fixed pitch angle k, the absorption coefficient at frequency 
v' < 0.37 2 (sinK)fg e = 0.37 2 (sin/t)(eB'/27rm ( ,c) is 
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in the bulk frame (e.g. Rybicki & Lightman 1979). 
tical depth t u = a' u {r/Y c ) oc j e , and so there is 
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The op- 
a critical 



energy %(t v = 1) above which the medium is optically thin 
at any given frequency. We normalize this frequency 12 to 
v 1 = T- l (l eV/h) ~ 5 x lO^iySO) -1 Hz in the bulk frame. 
Relating the pair density n'± to the Thomson optical depth t'±, 
one finds 
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Setting B' = B/T c = 2 x \tf B s (T c /50Y l G in the gamma-ray 
emission zone (we expect the magnetic field to have a weak 
dependence on L m \ andM w : eq. [74]) gives 
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In practice, the optical depth in such energetic particles is lim- 
ited to t± ~ 1 / 7 2 , and so we have 



7e ( r „ = 1) = 1 x lO 2 ^ 11 (smn) 2/u 
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(191) 



Synchrotron radiation in the optical band requires the elec- 
trons to be more energetic than 0.37 2 (sinK)/z(r c t^ e ) > 1 eV or, 
equivalently, 

le > Jemin = 50B" 1/2 (shm)" 1 / 2 . (192) 

The optical synchrotron emission is self-absorbed if 7 em j n < 
7e < lA T v = 1), and is determined only by the temperature 
k B Tj ~ |7 e m e c 2 of the electrons. The energy density of optical 
photons in the bulk frame is 
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The isotropic luminosity of the optical photons Qtv 1 = 1 eV/h) 
is 
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Here we have substituted expressions (38) and (39) for the de- 
celeration radius /?decei+ and the Lorentz factor at this radius. 

One possible way of avoiding excess optical/IR emission is 
to inject electrons with random Lorentz factors -f e ~ 10 4 - 10 5 , 
so that their synchrotron radiation is in the gamma-ray band 
(e.g. Rees & Meszaros 1994). The synchrotron output in the 
optical/IR band is then suppressed by a factor (/jw/£ pea k) 1//2 ~ 
10" 2 - 10" 2 5 . However, such a large mean energy per particle 
is inconsistent with the level of pair creation that must neces- 
sarily occur if decelerating agent is a dense Wolf-Rayet wind or 
breakout shell (§7.2). In such a compact environment, the mean 
energy per electron/positron is reduced to o± ^£' b /t± <~3-10. 
The cooling spectrum emitted by the particles (dL/dlnE^ oc 

1/2 

E 7 ' ) is also inconsistent with the observed low-energy spectra 
of GRBs (Ghisellini & Celotti 1999). 

The optical synchrotron emission is also suppressed if the 
high-energy portion of the gamma-ray spectrum arises from 
inverse-Compton scattering of thermal seed photons, and if 
shock acceleration is not the principal energization mecha- 
nism. The inverse-Compton spectrum extends up to an energy 
<~ 47 2 fcB?bb, which means that a distribution of particle energies 
up to 7 C < 10-30 is sufficient to create an extended gamma- 
ray spectral tail. The maximum energy to which electrons and 
positrons may be flash heated by decaying MHD turbulence is 
j e ^ o±, which is typically less than <~ 10 2 in the prompt decel- 
eration phase. Their synchrotron emission is also suppressed by 
a factor <~ 1 /7 2 In T c compared with the inverse-Compton emis- 
sion (in the case of a flat particle energy spectrum; eq. [159]). 

By contrast, the particle spectrum that is generated by shock 
acceleration is cut off by synchrotron cooling at a much higher 
energy, where the synchrotron energy is <~ d~^jn ( ,c 2 . First- 
order fermi acceleration is therefore disfavored as the mecha- 
nism for heating particles behind the reverse shock, given that 
a GRB radiates only < 10~ 3 of its bolometric output in the 
optical-IR band. The same conclusion applies to other accel- 
eration mechanisms which operate on a short timescale, includ- 
ing the resonant absorption of ion cyclotron waves by positrons 
(Hoshinoetal. 1992). 



8.5. Transition of Prompt Gamma-Ray Burst to Afterglow 

In the model advanced in this paper, the prompt gamma-ray 
emission continues as long as the shell of seed thermal photons 
overlaps the shell of turbulent magnetofluid. We now examine 
the transition from the prompt burst emission to the afterglow. 
A good example is GRB 980923 (Giblin et al. 1999). At such 
a transition, the prompt emission is often much more variable 
than the afterglow. 

Several factors point to the region between the reverse shock 
and the contact as the source of the prompt gamma-ray emis- 
sion. First, strong variability in the burst emission requires 
beaming of the radiation in the bulk frame (§7). Second, the 
radiative efficiency is larger when the pairs contribute a larger 
fraction of the particle inertia, which is the case behind the con- 
tact (compare eqs. [54] and [177]) And, third, the magnetic 
field is expected to be stronger between the reverse shock and 
the contact (§5.2). 

The Wolf-Rayet wind remains pair-loaded and mildly rela- 
tivistic beyond the point at which the reverse shock has com- 
pleted its passage through the ejecta shell. We have found that 

'2/3 —2/3 

it is still being pushed to r am b ~ 2M^_ 5 A/j at the radius 
^recei+ (eqs. [38] and [40]), independent of L re i. Beyond this 
point, r am b decreases inversely with radius. One infers that 
the ambient medium becomes fully quiescent only after the 
seed thermal photons have decoupled from the forward shock. 
The time since the burst scales as r/T 2 oc r 5 / 2 , and so one ex- 
pects an observable transition in the brightness of the X-ray 
synchrotron emission from the forward shock around a time 
t <~ 30(l+z)M w _5 s following a burst. Beyond this point, a 
much larger fraction of the energy density behind the forward 
shock is in slowly cooling ions than during the prompt burst 
phase. It should be emphasized that there is no clear connection 
between the amplitude of the prompt gamma-ray emission and 
the immediate post-flare synchrotron emission in this model. 

A rapid decline in the X-ray flux has been observed following 
several Swift bursts, the most plausible explanation for which 
is off-axis emission from the prompt emission phase (Taglia- 
ferri et al. 2005). Since the overlap of the prompt MeV pho- 
tons with the forward shock will not cut of sharply (we expect 
r T ~ 1 at the end of the prompt emission phase), this would re- 
quire strong dissipation of the non-radial magnetic field behind 
the reverse shock, if the ejecta shell were spherical. 

The contact discontinuity is, however, expected to be 
strongly nonspherical (§7.2). We have argued that Rayleigh- 
Taylor modes are excited in the cooling breakout shell over a 
wide range of angular scales. Relativistic ejecta penetrating 
through holes in the breakout shell will form 'subjets' that are 
boosted with respect to the shell material. The size of the holes 
decreases with the radius at which they are created, and the final 
generation of holes is created on an angular scale < T" 1 . One 
may therefore expect a few subjets to be present within each 
causal patch of area (r/T c ) 2 (cf. Kumar and Piran 2000). 

We note that some bursts show very smooth pulses which 
decay exponentially. These events appear to be dimmer on av- 
erage, and are rare amongst the brightest bursts (Norris et al. 
2005). Some possibly involve a different emission geometry, 
such as the breakout of a mildly relativistic shock from the sur- 
face of the Wolf-Rayet star (Tan et al. 2001) that is powered 
by a buried jet. It should also be noted that the optical depth 
in pairs grows toward lower values of the isotropic luminosity 
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L re i and higher Wolf-Rayet mass loss rates M w , because the rel- 
ativistic outflow is stopped more rapidly: t± cx L^/ 2 M„ 6 (eq. 
[177]). Variability in the burst light curve would be washed out 
when t± 3> 1 . 

There are suggestions from the analysis of broadband after- 
glow emission that the mass density in the ambient wind mate- 
rial is as low as M w -$ /V w $ ~ 10~ 3 - 10" 2 in a few bursts (Kumar 
& Panaitescu 2003; Chevalier et al. 2004). In this case, one can 
still define a radius at which the the photon shell will decouple 
from the relativistic ejecta (in the absence of a breakout shell). 
Taking the Lorentz factor of the ejecta to be r re i = Tb a it (eq. 
[12]) gives the decoupling radius 
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In this expression, we have made use of the relation Li so = 
^rei/(l _ £bb)- The compactness in the bulk frame is 
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Note that a much lower compactness would be obtained if the 
ejecta began to expand freely from a smaller radius ~ 10 7 - 10 8 
cm (rather than the photosphere R+ of the Wolf-Rayet star). It 
should also be emphasized that £' is sensitive to the asymptotic 
Lorentz factor of the ejecta shell, £' cx r~ 5 ,. This means that, 
if the reverse shock is sub-relativistic, then the photon shell 
will decouple from the ejecta shell before the reverse shock has 
completed its passage through the ejecta shell. However, a sub- 
relativistic reverse shock implies that r re i < SOL'/^, Afj" 1 ' 6 (eq. 

[39]), in which case £' > 6L re ]5 4 [ At] I/<6 at the radius of decou- 
pling. The portion of the magnetized shell which dissipates in 
the presence of the thermal photon bath will, therefore, do so at 
a large compactness. 

8.6. The Peculiar (?) Burst 941017 

The hard, high-energy spectral component that was observed 
in GRB 941017 (Gonzalez et al. 2003) has interesting implica- 
tions for the mechanism of particle heating. The energy index 
was observed to flatten from (3 ~ — | to (3 ~ at high energies. 
The high-energy component decayed more slowly, and so the 
transition energy between the two power-law components also 
decreased with time, from ~ 10 MeV down to ~ 300 keV. The 
burst had a very long duration, Too ~ 10 2 s, and its fluence was 
one of the dozen highest measured by BATSE. The source was, 
therefore, probably located at a cosmological redshift z < 1 . It 
is not clear how atypical this high-energy emission is, given 
the brightness of GRB 941017 and the relative difficulty of de- 
tecting it in lower-fluence bursts. The presence of strong high- 
energy emission in bursts with relatively soft 1-10 MeV spectra 
would significantly boost the pair creation rate and the optical 
depth in the emitting plasma. 

Inverse-Compton scattering of the ~ 0.5 MeV peak photons 
is the most obvious source of the high-energy component. If the 

13 Including the Klein-Nishina suppression of the scattering cross section. 



cooling particles are isotropically distributed, then the (3 = in- 
dex is inconsistent with passive Compton cooling ((3 = —j) and 
so the particles must be continuously heated (Stern & Poutanen 
2004). The slope of the rising spectrum would then mirror that 
of the spectrum below £p e ak; indeed (3 ~ is the most com- 
monly measured low-energy spectral slope in GRBs. 

This line of reasoning allows one to place some significant 
constraints on the mechanism of particle heating. In particu- 
lar, one can show that the heated particles must occupy a small 
fraction of the volume of the ejecta, during the peak of the burst 
emission. Given that the dominant contribution to the bolomet- 
ric output of GRB 941017 is at high energies, each electron or 
positron must, on average, radiate a total energy <~ a±m e c 2 in 
inverse-Compton photons. Note that the value of the magne- 
tization parameter a± in the high-energy emission zone may 
be significantly larger than the average value within the ejecta 
(§8.1.2). The cooling time at a random Lorentz factor ^ e > 30 
is much shorter than the burst duration, t coo \/At ~ i/£'j e ( as 
measured by the observer). This means that each particle 
must perform its emission duties over a very short timescale, 
femission/Af ~ (<7± / -f e )t cooi ~ <j±/£'j 2 < 10" 4 cr ± . Energy must, 
therefore, be transferred to the particles over a short, but not 
microscopic timescale. For example, shock acceleration does 
not satisfy this constraint, because the residency time of a par- 
ticle of energy j e <~ 30 near the shock is orders of magnitude 
shorter than f coo i. One infers that the outflow must have small- 
scale structure, on a scale < 10~ 4 er± cAt ~ (0.01 s)cr± c, which 
allows dissipation to take place so rapidly. Indeed, reconnec- 
tion of a reversing magnetic field in the outflow will have the 
desired effect (§5.1; 8.1.2). 

Could the (3 = spectrum be a signature of slow Compton 
cooling of the emitting particles? This requires the compact- 
ness to be tuned to a value somewhat smaller than ~ 7J 1 < 
0.03. During the peak of the burst emission, we deduce £' ~ 

-1/4 * 5/ 12 -5/12 

5L rel51 M w _ 5 V w8 from eq. (46), which indicates that the ex- 
ternal wind inertia must be exceedingly small for this explana- 
tion to hold. It is, nonetheless, clear that the compactness of 
the ~ 1 MeV photons must drop rapidly toward the end of the 
burst. The increasing dominance of the high-energy component 
is therefore consistent with a drop in r± and an increase in the 
cooling time (at fixed j e ) at the reverse shock. 

Synchrotron self-Compton emission by 7 c ~3x 10 4 elec- 
trons at the forward shock is another possible high-energy 
emission mechanism (Granot & Guetta 2003; Pe'er & Wax- 
man 2004a), which was motivated by the delayed onset of 
the high-energy component in GRB 941017 and the apparent 
need for slow cooling to explain the hard (3 = spectrum. We 
note that for the parameters e B < 10" 4 and T c ~ 300 adopted 
in this model, Compton cooling by the MeV photons will, 
in fact, dominate 13 cooling by optical synchrotron photons as 
long as the gamma-ray photons overlap the forward shock (Be- 
loborodov 2005). 

9. CONCLUSIONS 

We have examined, in some detail, the hypothesis that 
long gamma-ray bursts are emitted by relativistic material that 
breaks out from the photosphere of a Wolf-Rayet star, and 
then propagates through the surrounding dense wind. We have 
started with a key simplifying assumption, that there are only 
two major sites of dissipation in the outflow: where the rela- 
tivistic jet breaks out of the Wolf-Rayet star (at R t ~2x 10 10 
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cm); and the much larger zone where the reverse shock makes 
significant backward motion through the ejecta shell (at a dis- 
tance 10 14 - 10 15 cm). The key components of the relativis- 
tic outflow are a non-radial magnetic field and a thermal bath 
of hard X-ray photons. The magnetic field changes sign many 
times, due to the stochasticity of the dynamo operating in the 
central engine. The photons are created in the turbulent jet core 
as it is heated by Kelvin-Helmholtz instabilities during the pas- 
sage of the jet through the envelope of the Wolf-Rayet star. The 
kinetic energy of the entrained baryons can be comparable to 
the energy carried by these other two components, but does not 
have to be. 

Our analysis has covered the deceleration of the ejecta shell, 
as influenced by pair-creation in the ambient wind material; the 
entrainment of a thin shell of stellar material by the jet head; 
and the structure and composition of the fluid on either side 
of the contact discontinuity. We have described a mechanism 
by which turbulent motions in an expanding, relativistic mag- 
netofluid will transfer energy to electrons and positrons, and 
then (by Compton scattering) to an ambient radiation field. This 
inverse-Compton radiation is beamed along the magnetic field 
in the bulk frame, thereby providing an explanation for the rapid 
variability of many gamma-ray bursts, without resort to an ill- 
defined source of internal variability in the outflow. The ulti- 
mate source of seed photons is the relativistic jet itself (and not 
in its sub-relativistic cocoon) as it suffers Kelvin-Helmholtz in- 
stabilities during its emergence from the Wolf-Rayet star. 

The normalization and slope of the Amati et al. (2002) rela- 
tion between peak energy and isotropic energy of long gamma- 
ray bursts is reproduced, essentially without free parameters. 

1 li 

This relation is predicted to flatten from Speak ~ to E pei ± <~ 
Lj 1 ^ 4 below a luminosity L; so ~3x 10 5() ergs s" 1 . This transi- 
tion corresponds to a wide jet opening angle, > 1 /V3, at the 
surface of the Wolf-Rayet star. 

One also obtains an explanation for why the short gamma- 
ray burst GRB 050509b (Bloom et al. 2005) did not follow this 
relation: the relatively hard spectrum of this burst reflected a 
smaller radius at which the ejecta went into free expansion, and 
at which the temperature of the seed black body photons was 
frozen. The bright initial spikes of the SGR flares are another 
example of energetic bursts which deviate from the Amati et al. 
relation, and are consistent with a smaller thermalization radius 
of 10-100 km (Hurley et al. 2005; Palmer et al. 2005). 

We close with a summary list of some addition conclusions. 

9.1. Additional Conclusions 

1. The passage of the reverse shock wave through the 
ejecta shell occurs when the compactness of the radi- 
ation field near the contact discontinuity has dropped to 
a characteristic value £' ~ 10. At this point, the gamma- 
ray flux across the foward shock is no longer sufficient 
to push the ambient medium up to a speed comparable 
to that of the contact discontinuity. There is no detach- 
ment of the pair-loaded wind material from the ejecta 
shell. 

2. The prompt gamma-ray emission takes place in a 
medium with an optical depth ~ 1 in electron-positron 
pairs. Pairs generated by side-scattering of gamma-ray 
photons in the Wolf-Rayet wind are the dominant source 
of scattering particles in between the forward shock and 
the contact. Pairs are created at a somewhat higher rate 



by collisions between photons in the turbulent region 
between the reverse shock and the contact. This shocked 
material becomes optically thin only outside ~ 10 15 cm 
from the central engine. 

3. The shocked Wolf-Rayet wind and the shocked rela- 
tivistic outflow are themselves subject only to a weak 
Rayleigh-Taylor instability at their mutual contact dis- 
continuity. However, the shell of matter that is entrained 
by the jet head during the final breakout from the star 
is much denser. This breakout shell will gradually ac- 
celerate as it intercepts relativistic ejecta from behind, 
and will suffer from a strong corrugation instability as it 
cools. The radiation that is trapped by the breakout shell 
has a temperature of <~ 2 keV when the shell becomes 
optically thin. 

4. The interaction of the relativistic ejecta with the frag- 
menting breakout shell offers a hybrid mechanism for 
triggering dissipation and gamma-ray emission. In this 
case, the breakout shell does not form a truely 'external' 
medium, since its structure depends on its interaction 
with the relativistic wind material at a smaller radius. It 
is especially important when the mass-loss rate in the 
Wolf-Rayet wind preceding the burst is much less than 
10~ 5 M Q yr" 1 . The role of the corrugation instablity in 
generating long null regions in gamma-ray bursts should 
be investigated. 

5. Electrons and positrons in the outflow are electrostati- 
cally heated in regions of strong MHD turbulence. The 
Kolmogorov energy flux in torsional waves is damped at 
a smaller frequency than the ion gyrofrequency. We find 
that similar conclusion applies to the magnetic corona of 
a black hole, and to any magnetofluid in which B 2 /8tt is 
larger than the rest energy density in electrons (pairs). 

6. The radiation emitted by charges in a turbulent mag- 
netofluid is beamed along the background magnetic 
field. The damping of MHD turbulence in a medium 
with (B') 2 /Stt 3> n'±m e c 2 leads to electrostatic acceler- 
ation of charges along the magnetic field. Indeed, when 
the gamma-ray emission is triggered by the decelera- 
tion of the relativistic ejecta, one requires beaming in 
the bulk frame to explain the fast variability observed 
that is observed in many bursts. 

7. The outflow is photon rich (« 7 /n± ~ 10 2 - 10 3 ) and the 
photons are the dominant coolant of the electrons and 
positrons - even if the magnetic field has a higher energy 
density. Numerical simulations (e.g. Pe'er, Meszaros, & 
Rees 2005) which do not account for the anisotropy or 
spatial inhomogeneity of the electron distribution will 
lead to qualitatively different results. 

8. The observed correlation St cx £ 7 between photon en- 
ergy and the width of gamma-ray pulses arises natu- 
rally from inverse-Compton cooling of pairs that have 
been flash-heated up to a limiting Lorentz factor (7,,) <~ 
£' b (SB'/B') 2 /t± > 1. Some substructure in GRBs does 
not show clear asymmetries between rise and decay; we 
ascribe it to variations in the direction of beaming as- 
sociated with wave excitations of the background mag- 
netic field. 
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9. The electrons/positrons that are heated by decaying 
MHD turbulence behind the reverse shock have a nega- 
tive gradient in energy away from the shock. There is, 
correspondingly, a gradient in the energy of the inverse- 
Compton photons. The harder spectra that are often ob- 
served at the beginning of GRB sub-pulses could be ex- 
plained by this effect. 

10. Synchrotron emission in the optical range is suppressed 
because the heated electrons and positrons have small 
pitch angles. The inverse-Compton emission of MeV- 
GeV gamma-rays only requires particle energies up 
to 7e ^ 30- 10 2 , which is smaller than the energy 
that is needed for optical synchrotron emission Qy e <~ 
50(sinK)" 1//2 ) when the small pitch angle is taken into 
account. 

11. The magnetization is very different in the fluids that 
have passed through the foward and reverse shocks: 
the magnetic field probably contributes a much smaller 
fraction of the pressure in front of the contact. The 
Weibel instability generates a very small-scale magnetic 
field. This field has a tendency to smooth out and in- 
crease its scale downstream of the shock (Medvedev et 
al. 2005). However, if the rate of smoothing is at least 
<~ 10" 4 of the maximal rate allowed by causality, then 
the r.m.s. smoothed field will be weaker than the field 
which is swept up from the Wolf-Rayet wind and lin- 
early compressed behind the forward shock. 

12. The thermal photons advected out from the sub- 
photospheric region of the Wolf-Rayet star have a finite 
duration, and a transition from prompt burst to afterglow 
can be explained by a drop in the thermal photon den- 
sity. A modest anisotropy in the gamma-ray emission 
on a angular scale <~ T~ l , due to corrugation instabili- 
ties of the breakout shell, could also influence the shape 
of the transition to the afterglow regime. 

13. The gamma-rays themselves are the dominant coolant 
at the foward shock. As pre-acceleration turns off, the 
medium ahead of the shock develops relativistic mo- 
tion with respect to the contact discontinuity (a differ- 
ential Lorentz factor <~ 20). The pairs which pass across 
the forward shock emit inverse-Compton gamma-rays 
with an energy ^T± <~ 500 times the peak energy of the 
burst. This emission is, however, too spectrally soft and 
too weak to explain the rising high-energy component 
of GRB 941017. We argue that the high-energy emis- 
sion in GRB 941017 requires rapid (but not too rapid) 
heating of the particles behind the reverse shock, most 
probably triggered by magnetic reconnection. 

14. The global MHD instabilities identified by Lyutikov and 
Blandford (2003) are effective at tangling the magnetic 
field while the jet is still working its way through the 
envelope of the Wolf-Rayet star. They are probably less 
important in the prompt emission zone at ~ 10 14 - 10 15 
cm. The Lorentz factor of the contact is still large at 
this stage, T c <~ 50, and much larger than the inverse 
of the jet opening angle. Smaller angular structures 
in the outflow will be created by the corrugation in- 
stability of the breakout shell; but the growing mode 



will have a fairly large angular scale close to the star, 
9 ~ r~' oc (r//?*)" 1 / 3 , where the Lorentz factor of the 
contact is still modest. The angular scale decreases as 
the breakout shell is accelerated outward, but at the end 
of the prompt deceleration phase the contact is moving 
nearly as fast as the ejecta. 

We close this paper with a question: What is the main out- 
standing theoretical problem relating to the prompt emission 
of gamma-ray bursts? This question is far from being rhetor- 
ical: the proper design of numerical experiments to study the 
3-dimensional behavior of plasmas - especially in the extreme 
regimes that are encountered in GRB outflows - requires a basic 
understanding of the dominant physical effects that are likely to 
be encountered. 

Most recent research would seem to point to the following 
answer: one must understand the structure and intermittency of 
the outflow from a stellar-mass black hole that is surrounded by 
a neutronized torus. However, this problem is vast, as it encom- 
passes the mass flow in the torus; the Blandford-Znajek process 
operating in the black hole magnetosphere; the acceleration and 
collimation of the outflow near the engine; the interaction of 
this jet with the envelope of the Wolf-Rayet star; nuclear pro- 
cesses occurring over a wide range of radii; and the various 
instabilities that occur in the jet outside the star. In this picture, 
the emission problem cannot simply be reduced to the problem 
of understanding particle acceleration and radiation at a mildly 
relativistic shock wave, since the details of the microphysics 
depend strongly on where the shock occurs in the outflow, the 
strength of the magnetic field, and the feedback of radiation 
processes on the particle composition and energy distribution. 

We are suggesting here that, in fact, the theoretical problem 
is a bit more tractable. Because the non-thermal emission is 
triggered by the interaction with the external wind and the ma- 
terial of the breakout shell, one can focus down on a small part 
of the outflow, and in particular on the problem of how long- 
wavelength magnetohydrodynamic turbulence is damped in a 
fluid satisfying the following conditions: approximate equipar- 
tition between thermal radiation and magnetic field; B 2 /8tt > 
n±m e c 2 ; and compactness I <~ 10. Because the dissipation zone 
has been localized to a well-defined region, one can now begin 
to examine the radial structure of the outflow in a constrained 
way, and to understand the origin of distinct emission compo- 
nents such as was observed in GRB 941017 (Gonzalez et al. 
2003). 

From this perspective, the classical gamma-ray bursts have 
a greater affinity with the giant flares of the Soft Gamma Re- 
peaters that has been generally assumed in the recent literature. 
Much of the early confusion about whether gamma-ray bursts 
originate in the magnetospheres of Galactic neutron stars, or 
in relativistic fireballs at cosmological distances, turns out to 
be connected to a similarity in the mechanism underlying these 
spectacular releases of radiant energy. 
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